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ABSTRACT 


This  report  is  the  first  of  a  two-year  effort  for  deter¬ 
mining  the  dynamic  and  static  aerodynamic  stability  derivatives 
for  volute  stabilized  cylindrical  forebodies.  The  complete 
analytical  treatment  of  both  rigid  and  flexible  models  is  dealt 
with  and  these  results  compared  to  data  obtained  from  the  Eglin 
low  speed  wind  tunnel.  The  theory  presents  several  methods  that 
give  the  designer  necessary  techniques  for  estimating  Cma  and 
(Cinq  +  Cm&)  for  a  variety  of  different  volute  shapes  and  cylin¬ 
drical  forebodies.  These  techniques  yield  reasonable  answers 
for  both  rigid  and  flexible  models  operating  at  Rn  =  2  x  105 
and  Moo  =  0.2  with  volute  tails  sufficiently  long  so  that  the 
potential  flow  theory  used  is  valid.  The  complete  data  obtained 
at  both  the  Eglin  facility  and  at  Arnold  Engineering  Development 
Center,  Arnold  Air  Force  Station,  Tennessee,  joined  with  the 
analysis  presented,  served  to  show  that  the  volute  provides  an 
effective  yet  compact  stabilizer  for  Rn  =  2  x  105  to  1  x  10° 
for  Moo  =  0.2  to  0.5  for  cylindrical  forebodies. 


Distribution  limited  to  U.  S.  Government  agencies  only;  this 
report  documents  tests  and  evaluation  of  potential  military 
hardware;  distribution  limitation  applied  May  1971.  Other 
requests  for  this  document  must  be  referred  to  the  Air  Force 
Armament  Laboratory  (DLRA) ,  Eglin  Air  Force  Base,  Florida 
32542. 


iii 

(The  reverse  of  this  page  is  blank) 


TABLE  OF  CONTENTS 


Section  Title  Page 

I  INTRODUCTION .  1 

II  RIGID  BODY  ANALYSIS .  4 

HI  FLEXIBLE  BODY  ANALYSIS .  37 

IV  FINAL  DESIGN .  68 

V  CONCLUSION .  73 

Appendix 

I  ANALYSIS .  75 

II  DATA .  97 


REFERENCES 


123 


s 


s' 

f 

i 

i 


i 


i 


LIST  OF  FIGURES 


Figure  Title  Page 

1  Volute  Stabilizer  .  1 

2  lv  Versus  xt  @  0  *  0 .  8 

3  lv  Versus  xt  @  9  =  0 .  9 

4  lv  Versus  ht  @  0  =  0 .  10 

5  lv  Versus  Xyt  @  0  =  0 .  12 

6  6  Versus  Sy/cT2 .  15 

7  0  Versus  Cms  Model  3+1 .  18 

8  6  Versus  Cms  Model  3+2 .  19 

9  0  Versus  Cms  Model  3+3 .  20 

10  0  Versus  Cms  Model  3+5 .  21 

11  0  Versus  Cms  Model  3+6 .  22 

12  0  Versus  Cm(X  Method  2  . .  23 

13  0  Versus  CmQ[  Method  3  .  24 

14  Sy/St  Versus  Cnto  .  27 

15  Master  Design  Curve .  29 

16  ly  Versus  0triin  . . 32 

17  ly/dy  Versus  (Cmq  +  Cm£)  .  34 

18  lv  Versus  (Cmq  +  Cm£) .  36 

19  lv  Versus  Cmaf .  43 

20  6  Versus  Cms  Flexible  Model  1 .  45 

21  0  Versus  Cms  Flexible  Model  2 .  46 

22  d2  Versus  di  .  60 

23  lv  Versus  6  .  63 

24  lv  Versus  8  .  64 

25  lv  Versus  Frequency  (Short  Period) .  65 

26  lv  Versus  Frequency  (Long  Period) .  66 

27  lv  Versus  Cmaf  . 72 


v 


vi 


[ 


LIST  OF  FIGURES  (Concluded) 


Figure  Title  Page 

1 1  - 1  Rigid  Volutes:  1,  2,  3 .  99 

1 1  - 2  Rigid  Volutes:  4,  5,  6 .  100 

1 1  -  3  Cylinder  Models  .  101 

1 1  - 4  Helical  Coil  Flexible  Volutes  .  102 

1 1  -  5  0  Ver  ses  Cn  For  Rigid  Volutes  .  103 

1 1  -  6  0  Verses  C^  For  Rigid  Volutes  .  104 

1 1  -  7  0  Verses  Cms  For  Rigid  Volutes  .  105 

1 1  -  8  0  Verses  Cn  For  Rigid  Models  .  106 

1 1  -  9  0  Verses  C^  For  Rigid  Models  .  107 

11-10  0  Verses  C^For  Rigid  Models  .  108 

1 1  - 1 1  0  Verses  C^  For  Flexible  Models  .  109 

11-12  0  Verses  C^For  Flexible  Models  .  110 

11-13  Accelerometer  Stand  .  Ill 

11-14  Accelerometer  Output  Filter  .  112 

11-15  Dynamic  Shaker  Stand  .  113 

11-16  Dynamic  Shaker  Stand  (Top  View)  .  114 

11-17  Von  Karman  Data  Sample.  0.75  Inch  Cylinder 

With  Small  Volute  . .  115 

11-18  Von  Karman  Data  Sample,  0.75  Inch  Cylinder 

With  Small  Volute .  116 

11-19  Von  Karman  Data  Sample,  1.25  Inch  Cylinder 

With  Small  Volute  .  117 

11-20  Von  Karman  Data  Sample,  1.25  Inch  Cylinder  With 

Small  Volute  .  118 

11-21  Volute  Tails  USed  In  Von  Karman  Test  Program  119 


LIST  OF  TABLES 


Table 


Title 


Page 


I  Cnia/Rcicl . . .  17 

II  ht/l0MATRIX  .  26 

III  (Cmq  +  Cm<j)/Rad.  (THEORETICAL)  .  33 

IV  (Cmq  +  Cm^)/Rad .  35 

V  Cma^/Rad . 40 

VI  BEST  DESIGN  LENGTH  VOLUTE  .  42 

VII  Cmaf/Rad .  44 


VIII  TRIM  ANGLE  . 

IX  (Cmq  +  Cm<.)  /Rad  (THEORETICAL) 

X  Aij  . 

XI  Cma/Rad . 

1 1  - 1  Volute  Model  Measurements  .... 

1 1 - 2  Cylinder  Model  Characteristics 

1 1  -  3  Volute  Stabilization  Data  ... 


47 

48 
56 
70 
98 

119 

120 


viii 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 


hfc  -  Moment  arm  measured  forward  from  the  model  eg  to  the 

center  of  pressure  of  cylindrical  portion  forward  of  the 
model  eg. 

• 

hb  -  Moment  arm  measured  aft  of  the  model  eg  to  ne  center  of 
pressure  of  the  cylindrical  portion  of  the  model  aft  of 
the  eg. 

h^  -  Moment  arm  measured  from  the  model  eg  to  the  center  of 
pressure  of  the  volute  stabilizer. 


ht  -  Moment  arm  measured  aft  of  the  model  eg  to  the  center  of 
pressure  of  aft  portion  of  the  model. 

hv  -  Moment  arm  measured  from  the  volute  hinge  line  to  the 
volute  center  of  pressure. 

A 

xt  -  Accounts  for  the  rotary  velocity  field  aft  of  the  eg. 


*b 


xb 

xVt 


Sb 

Sb 


-  Accounts  for  the  rotary  velocity  field  forward  of  the 
model  eg. 

-  Accounts  for  the  rotary  velocity  field  for  the  aft  por¬ 
tion  of  the  cylinder  relative  to  the  model  eg. 

-  Accounts  for  the  rotary  velocity  about  the  volute  mea¬ 
sured  from  the  model  eg. 

-  Accounts  for' the  rotary  velocity  field  about  the  volute 
measured  from  the  volute  hinge  line. 

-  The  planform  area  of  the  cylindrical  portiqn  of  the  model 
measured  forward  of  the  model  eg. 

-  The  planform  area  of  the  cylindrical  portion  of  the  model 
measured  aft  of  the  model  eg. 


Sv  -  Planform  area  of  the  volute. 
st  *  st  =  sb  +  Sv 

S  -  Total  planform  area  of  the  model, 

lb  -  Forebody  length 

lv  -  Volute  length 


ix 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS  (Continued) 


10  -  Length  of  the  cylinder  aft  of  the  model  eg. 

xCg  -  Model  center  of  gravity  measured  relative  to 

model  nose. 

x  -  Variable  running  along  the  model  longitudinal 

axis . 

dc  -  Diameter  of  the  cylinder. 

dv  -  Diameter  of  the  volute. 

dy  -  Diameter  of  the  truncated  portion  of  the 

volute. 


q 

cnt0 

Cnb0 

Cnv0 

Cm 

Cms 

cma 

(cmq  +  CII1(j) 
cm<s 


Dynamic  pressure. 

Normal  force  coefficient  of  the  area  aft  of 
the  model  eg  at  90  degrees. 

Normal  force  coefficient  of  the  cylindrical 
portion  of  the  model. 

Normal  force  coefficient  of  the  volute. 

Total  pitching  coefficient. 

Static  portion  of  Cm. 

Pitching  moment  slope  coefficient. 

Damping  coefficient. 

Represents  the  static  pitching  contribution  of 
the  volute. 


Cm6 

Cmaf 

(Cmq  +  Cm£)f 
J 

A 

0 


Represents  the  amount  of  intrinsic  damping 
power  of  the  volute. 

Flexible  pitching  slope  coefficient. 

Flexible  damping  coefficient. 

A  number  which  proportions  the  effects  of 
flexible  volutes  as  related  to  rigid  ones. 

Proportionality  constant  for  flexible  analysis. 

Proportionality  constant  for  rigid  analysis. 


x 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS  (Continued) 


AB  -  Forebody  contribution  to  the  Cmct  equations. 

AB'  -  Forebody  contribution  to  the  (Cm<^  +  Cm^)  equation. 

V  -  Represents  the  moment  developed  by  the  volute. 

kv  -  Lateral  spring  constant  of  the  volute. 

aij  -  Coefficients  of  the  two-degree  angular  freedom  equa¬ 
tions. 

Ib  -  Mass  moment  of  inertia  of  the  cylinder  about  the  model 
pitch  axis . 

I  -  Total  body  moment  of  inertia 

Kv  -  Mass  moment  of  inertia  term  resulting  from  two-degree 
angular  freedom  moment  summation. 

Iv  -  Mass  moment  of  inertia  of  the  volute  measured  about  the 
hinge  line. 

Mv  -  Mass  of  the  volute. 

Rv  -  Distance  from  the  pitch  axis  to  the  hinge  axis. 

rv  -  Distance  from  the  volute  hinge  line  to  the  volute  mass 

center . 

wv  ■  Volute  weight. 

Voo  "  Relative  velocity. 

Cij  -  Coefficients  on  the  characteristic  equation  for  the 
flexible  models. 

0  -  Angle  subtended  by  the  free  stream  velocity  vector  and 

the  model  longitudinal  axis. 

p  -  Air  density. 

s  -  Characteristic  equation  variable. 

dj_  -  Coefficients  of  the  quadratic  factors  of  the  4th  degree 
characteristic  equations. 

Pij  -  Elements  in  the  Routh-Hurwitz  stability  scheme. 


xi 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS  (Concluded) 


R  -  Cylinder  cavity  radius, 

t!  -  Volute  spring  material  thickness. 

t2  -  Volute  spring  material  width, 
n  -  Number  of  volute  helical  coils. 

60  -  Volute  half-angle. 

6  -  Deflection  of  volute  relative  to  cylinder  longitudinal 

axis  . 


SECTION  I 


INTRODUCTION 


This  report  documents  the  results  of  investigations  of 
the  mechanical  behavior  of  volute  stabilized  cylindrical  fore¬ 
bodies.  The  cylinder  model  flight  characteristics  were  deter¬ 
mined  through  experimental  and  analytical  programs  which  yielded 
values  for  the  stability  derivatives  and  several  predictive 
schemes  for  evaluating  Cma,  (Cmq  +  Cmd )  • 

The  volute  stabilizer  is  essentially  a  tapered  helical 
spring  which  can  be  packaged  conveniently  in  a  small  area  at 
the  base  of  a  forebody  (Figure  1).  Spring  steel,  which  provides 
strength  and  flexibility,  is  an  ideal  material  for  volute  stabi¬ 
lizer  construction. 


Figure  1.  Volute  Stabilizer 
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The  basic  dimensionless  ratios  are:  ly/dy,  dy/dc,  lv/ l** 
TJie  volute  is  considered  to  be  fully  characterized  once  ly,  b 
dv,  dv,  ky  (spring  constant)  wy,  and  ly  have  been  specified. 

For  a  given  volute  planform  area  and  forebody,  the  volute 
stabilizer  is  considered  effective  if  Cma  and  (Cmq  +  Cma)are 
negative.  The  determination  of  these  derivatives  as  functions 
of  forebody-volute  dimensions  represents  the  core  of  the  fol¬ 
lowing  sections.  Once  they  are  determined,  the  frequency  equa¬ 
tion  can  be  solved  for  the  characteristic  roots  of  the  model. 
The  determination  of  Cma  and  (Cmq  +  Cma)is  complicated  by  the 
fact  that  volutes  are  flexible  and  are  forced  to  operate  in  the 
slipstream  of  the  forebody. 

Because  of  the  success  of  linear  aerodynamics  in  predict¬ 
ing  stability  and  motion  for  small  coordinate  perturbations  of 
many  shapes,  the  first  approach  to  the  problem  was  to  force  the 
nonlinear  rigid  body  equation  of  motion  Equation  (5)  into  a 
linear  form.  The  basic  moment  summation  required: 


Constrained  eg,  pitch  only,  18  =  qoodSCm  (1) 

.  d 

Where,  Cm  =  Ctrim  +  Cma  9  +  (cmq  +  cmd)02V^  (2) 

Cm  was  derived  in  terms  of  dimensions  of  the  forebody-volute 
combination  and  the  normal  force  coefficient  Cn.  Nonlinear 
effects  were  later  entered  into  the  basic  form  by  allowing  for 
forebody  interference  and  sin20  variation  of  the  normal  force. 
The  results  of  the  rigid  body  analysis  are  divided  into  three 
parts  corresponding  to  three  approaches  used  to  solve  for  Cm. 
The  first  is  a  linear  theory,  the  second  includes  forebody 
interference,  the  third  involves  interference  and  sin20  varia¬ 
tion  of  Cm. 

Because  the  volute  tail  is  flexible,  serious  questions 
arise  as  to  whether  a  volute  equipped  cylinder  will  be  dynami¬ 
cally  stable.  Insuring  dynamic  stability  is  a  much  more  pro¬ 
found  problem  than  the  simpler  problem  of  guaranteeing  static 
stability,  criteria  for  which  are  set  up  in  the  rigid  body 
analysis.  Dynamic  stability  may  be  inferred  from  the  following 
set  of  equations  which  describe  an  elastic  model  having  only 
two  degrees  of  angular  freedom: 


eg  moment,  at  l  0+ai  2$+a1 30+a1 40+ai  5  |5+a1 6S  =  0  (3) 

Hinge  moment,  a2 19+a22^+a23Q+a2^+a2  s^+a2  66  =  0  (4) 
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The  entire  flexible  section  deals  with  these  equations.  It  is 
shown  in  the  analysis  that,  depending  upon  the  mass  and  inertia 
of  the  elastic  tail,  characteristic  equations  (frequency  equa¬ 
tions)  may  be  second  or  fourth  order  polynominals .  Stability 
is  investigated  by  invoking  the  Routh-Hurwitz  criteria. 

Basic  supporting  data  was  obtained  from  both  static  and 
dynamic  test  stands.  Dynamic  tests  utilize  two  special  stands 
developed  under  the  contract:  The  first  utilized  a  piezoelec¬ 
tric-equipped  model  to  drive  an  oscillograph  which  gave  accel¬ 
eration  versus  time  data  and  the  second  forced  the  test  model 
into  resonance  (forced  oscillation  testing).  Free-flight  data 
was  also  gathered  and  reduced. 
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SECTION  II 


RIGID  BODY  ANALYSIS 


The  rigid  model  equation  of  motion  is  obtained  from  a 
simple  moment  summation  about  the  forebody-volute  mass  center. 
As  shown  in  Appendix  I,  the  relative  velocity  is  the  vector 
sum  of  the  free  stream  velocity  and  a  rotary  velocity  field 
which  coexists  with  any  rotational  motion.  The  entire  stabi¬ 
lizing  aerodynamic  moment  is  equated  to  the  normal  pressure 
force  multiplied  by  an  effective  moment  arm.  The  resulting 
equation  is  nonlinear: 


"  10  =  ^pCn.£^h-j-  +  xt®)|Voo®  +  xt®  1 

"  Hpcnb0hb  (v“e  '  xb^lV^e  "  xb§|Sb 


(5) 


It  is  shown  that  Equation  (5)  reduces  to  a  linear  form  if  a 
rumber  of  assumptions  are  made.  When  this  is  done: 


-IS  *  qdc^ 


2cnb  *thtst 


V. 


dc2S 


l  + 


Cnt0^tst 


dcS 


cnbQ^b^bsb 

cnto^t^tst 


/\  dp  • 
e - SL  0 

2Voo 


Cnb0xbSb 

cntnxtst 


60 


Linear  aerodynamics  require: 


(6) 


ie 


0  +  Cma6 


where  , 


Cmtrim  0 


(7) 


By  identification  of  the  common  terms  in  Equations  (6)  and  (7) 


CnVo  and  Cnb0  are  the  normal  force  coefficients  of  the  volute 

and  cylinder  evaluated  at  90°.  Cnvo  decreases  with  an  increase 
of  the  ratio  lv/dv  (Appendix  II).  The  final  results  have 
justified  approximating  Cny0  and  t’nbo  average  values  for 

volutes  having  lv/dv  ranging  from  0.78  to  4.79  and  cylinders 
with  lc/dc  ranging  from  1.48  to  2.36.  Cnv  and  CnbD  will  be 
replaced  by: 


The  remaining  variables  are: 

Sb  -  Planform  area  of  the  cylinder  forward  of  the  eg. 
Sb'-  Planform  area  of  the  cylinder  aft  of  the  eg. 

Sv  -  Planform  area  of  the  volute. 

St  -  Equal  to  S^>  * +  Sv  . 

S  -  Equal  to  the  total  model  planform  area. 
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Cn^o  -  Equal  to  0.833. 

Cnt  -  The  effective  normal  force  coefficient  of  the  tail 
area  aft  of  the  eg.  Sv 

Cnt0  =  0.267  +  0.833. 

ht  -  Moment  arm  used  to  account  for  the  adverse  moment 
set  up  by  that  portion  of  the  cylinder  forward  of 
the  eg.  hb  =  *sxcg 

ht  -  Moment  arm  used  to  account  for  the  stabilizing 

moment  set  up  by  that  portion  of  the  cylinder  aft 
of  the  eg  plus  the  entire  volute,  ht  is: 


Xb  -  Linear  dimension  used  to  proportion  the  magnitude 
of  the  rotary  velocity  field.  It  accounts  for 
the  field  forward  of  the  cylinder  eg  giving  an 
approximate  value  to  (x)0  .  xb  =  ^cg* 

xt  Serves  the  same  purpose  as  Xb.  xt  is  complicated 
because  it  deals  with  two  different  planforms, 
i.e.,  rectangular  and  triangular.  For  all  analyti¬ 
cal  work  to  follow  it  is  defined  to  be  xt  *  -jd-v  + 


The  evaluation  of  ht  and  xt  is  actually  more  involved  than 
is  indicated  by  their  simple  definitions.  The  assigned  values 
followed  a  detailed  analysis  intg  the  behavior  of  each  with 
variations  of  lv/lo>  e>  and  0-  xt  was  found  by  evaluating  the 
following  integral  equated  to  an  average  expression: 


The  integration  gives #a  complex  equation  revealing  the  depen¬ 
dence  of  xt  on  9  and  fi  .  The  primary  interest,  however,  is  in 
the  range  over  which  $t  can  vary.  This  was  found  by  setting 
first  0  =  0  and  theij  6  =  0  in  the  results  of  Equation  (14). 
When  this  was  done  xt  was  found  to  obey: 


xt 


(15) 


6 


2  - 


TIq  3 


Tlpl2w+  lylo2  +  Tly3  0 


=  0 


210  +  lv 


(16) 


These  equations  represent  the  upper  and  lower  bounds  over 
which  xt  may  vary.  The  value  for  xt,  defined  earlier  and  used 
in  deriving  other  expressions,  differs  from  Equation  (15)  or 
Equation  (16)  (i.e.,  ido  +  — lv) .  xt ,  as  defined,  represents  a 

compromise  between  the  two  extremes  at  6  =  0  and  9=0.  How¬ 
ever,  it  may  be  necessary  for  some  design  problems  to  have  a 
better  value  for  xt,  and  for  this  purpose  Equation  (15)  and 
Equation  (16)  are  graphed  (Figures  2  and  3).  For  a  given  10, 
xtmay  be  found  as  the  volute  changes  length  (lv  varies)  for 
0=0  and  0=0.  The  same  argument  shows  that  ht  likewise 
varies  with  0  and  0.  The  integral  to  be  evaluated  does  not 
follow  directly  because  of  the  way  in  which  ht  was  originally 
derived,  ht  was  found  after  a  moment  summation  about  the  model 
eg  which  took  into  account  only  the  normal  force  acting  aft  of 
the  eg.  The  equation  is: 

CntQ^tSt  =  ^nb^b  Sb  +  Cnv^ht  Sv  (17) 

Eyidently  ht  is  a  function  of  two  other  moment  arms  hb  and  ht. 
hb  offers jno  problem  since  to  good  accuracy  it  can  be  approxi¬ 
mated  by  y10. 


o 


(18) 


Hqwever,  h^  does  vary  with  0  and  0.  ht  varies  with  ht  and 
ht  is  found  by  integrating  over  the  volute.  The  result  is: 

(19) 

(20) 


h'  -  1 


l0lv  +  7  + 


3  +  il3 


,0  =  0 


ht  =  1q  +  3  ^v ’ ®  =  0 


ht  is  graphed  (Figure  4)  for  values  of  lo  and  lv.  Equation  (17) 
with  the  corrected  ht  should  be  used  in  lieu  of  Equation  (13) , 
if  the  extra  accuracy  is  deemed  necessary.  Some  guidelines  on 
this  point  will  be  provided  in  the  discussion  on  stability. 
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Figure  3.  lv  Versus  xt  2  0  =  0 
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Although  it  is  not  used  in  the  linear  expressions  presented 
thus  far,  Xvt  serves  a  similar  purpose  in  future  work  as  xt.  It 
deals  with  tne  rotary  velocity  field  over  the  volute  only 
(whereas  xtAdealt  with  the  whole  tail  portion  aft  of  the  eg) . 
Values  for  xvt  are  (Figure  S)  ; 


x  yx  =  1q2  *  loly  *  lv2/ 2  0  =  0 

L  3  5 


*vt  =  xo  +  3xv,e  =  0 


Detailed  work  concerning  derivation  of  all  the  above  forms  for 
xt,  ht ,  and  xVt  are  included  in  Appendix  I. 


The  first  thing  that  is  apparent  from  Equation  (8)  is  a 
convenient  stability  criteria.  For  static  stability  it  is  re¬ 
quired  that  Cmct  <  0  or  that: 


1- 


CnbQhbsb  > 
CntohtSt 


0 


(23) 


which  is, 


htSt  >  ^  hbSb 
Lnt0 


The  largest  that  0.833/Cnto  can  be  is  1.00  so  that  Equation  (22) 
may  also  be  written  as: 


htSt  >  hbSb  (24) 


It  is  required  by  Equation  (24)  that  the  tail  planform  area 
and  ht  multiplied  must  be  greater  than  the  analogous  parameters 
hbSb  measured  forward  of  the  eg.  For  many  combinations  of 
cylinders  and  volutes  it  is  possible  to  have  St  =  Sb  in  which 
case  stability  depends  upon: 

ht  >  hb  (25) 

The  importance  of  ht  is  evident.  Normally  the  cylinder-volute 
combination  will  be  designed  to  insure  that  St  >Sb,  otherwise 
stability  may  be  compromised  by  variations  of  ht  with  e  or  0 
Large  St  relieves  the  necessity  of  having  precise  values  for  ht. 
For  determining  the  angular  frequency  and  damping  exponent 
(  e  -  3t  ) ,  use  the  best  obtainable  value  for  ht .  Equation  (17) 
will  serve  this  purpose. 
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Because  of  the  method  in  which  Equation  (5)  was  derived, 
it  is  impossible  for  the  rigid  model  to  be  dynamically 
unstable.  This  is  evident  from  Equation  (9)  which  can  never 
be  greater  than  zero.  TheAequation  is  nearly  the  same  as 
Equation  (8)  but  includes  x-t  and  xt).  Accurate  determination 
of  these  variables  will  be  manifest  in  3,  the  damping  constant. 
Dynamic  or  static  stability  of  the  rigid  model  does  not  depend 
on  Xf  As  will  be  shown,  this  is  not  the  case  for  the  flexible 
model . 

Before  discussing  the  numerical  results  of  the  preceding 
analysis,  two  additional  methods  for  computing  Cma  and  sub¬ 
sequently  Cms  will  be  presented.  Although  the  linear  analysis 
gives  good  answers,  it  fails  to  indicate  the  interference 
effects  from  the  forebody  and  nonlinear  variation  of  the  normal 
force  coefficient.  Returning  to  Equation  (5),  it  is  evident 
that  there  exists  a  82  term  which  approximates  sinze.  Writing 
the  exact  expression  for  the  static  pitching  moment: 


cms 


sine | sin0,-9O  <e<  90 


(26) 


The  pitching  moment  coefficient  in  linear  form  follows  directly 
from  Equation  (26)  and  is  again: 


The  pitching  moment  slope  coefficient  is: 


(27) 


3Cms 
3  0 


Cm(X  i .  e  .,  [Equation  (8)] 


(28) 


To  account  for  forebody  blanking  of  the  volute,  the  area  Sv 
must  be  altered  to  reflect  a  reduction  in  the  tail  surface 
which  can  act  to  stabilize  the  forebody.  The  detailed  deriva¬ 
tion  is  in  Appendix  I.  The  results  are: 

Sv  _ly 
d2  2d 

d  =  dv 


d 

81 


[tan0 


d  ]2 


(29) 


V 
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Figure  6  dnd  Equation  (29)  represent  the  variation  of  Sv  =  Sy 
with  0,  if  the  area  planform  is  a  triangle  (dy  -  0).  Most 
volutes  end  at  the  forebody  blunted  and,  therefore,  the  planform 
is  trapezoidal.  This  can  be  accounted  for  by  adding  to 
Equation  (29)  : 


where 


attack  is: 


dv'  (lv  -  xO 

(30) 

(dr  -  d„0 

tan6  +  dv  -  dv' 

(31) 

21y 

wind  stream  for 

a  given  angle  of 

-  q  t  +  q  H 
“  *>V 

(32) 

Sv  will  be  good  enough  to  approximate  the  area  if  dv  is  not 
large  or  not  more  than  about  1/10  dv.  Sy  is  plotted  for  dif¬ 
ferent  values  of  lv/cF.  To  use  the  new  area  in  the  Cm  expres¬ 
sion,  an  alternate  form  of  Equation  (8)  has  to  be  used  and  is: 


-m. 


a 


C 


ny0htSV 

-c'cir  ~~ 


+  cnbohbsb' 
CnvQ^t'Sy 


cnbQhbsb 

CnVoht£v 


(33) 


This  expression  for  C]tia  is  not  new  since  Equation  (8)  was  derived 
directly  from  it.  The  new  variables  are: 

/si  /s 

Xb  -  Serves  the  same  purpose  as  did  Xb  but  is  equal  to 
o  • 

hb  -  Moment  arm  for  the  normal  force  which  act^  on  that 
portion  of  the  cylinder  aft  of  the  eg.  hb  -  yl0 

Sb  -  Planform  area  of  the  cylinder  aft  of  the  model  eg. 

The  new  stability  criteria  is: 

(cnVQhtsv  +  cnb0hb'sb'>  >  CnbohbSb  (34) 

at  0  =  0  the  criteria  reduces  to: 

cnbQhbsb  >  Cnbchbsb  (35) 
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If  Equation  (35)  is  not  satisfied,  the  inference  is  that  the 
model  will  tend  towards  static  instability.  This  is  partially 
correct.  As  soon  as  the  model  tends  to  overturn,  Equation  (34) 
becomes  valid.  Since  Sv  rises  rapidly  with  0, increase  in  0 
immediately  tends  to  restore  stable  motion.  The  length  of  the 
volute  is  the  critical  parameter  since  short  volutes  do  not 
have  the  sharp  rise  in  Sy  exhibited  by  long  volutes  (Figure  6). 
With  the  new  values  for  Sv,  Cms  becomes  a  function  of  0.  With 
this  new  form  for  Cms  >  two  nonlinear  theories  are  evident. 
Equation  (26)  i§  the  most  nonlinear  form  investigated.  Replac¬ 
ing  sin20  with  0sin0  yields  a  third  form  which  is  a  compromise 
between  the  linear  theory  and  the  sin20  theory.  The  three 
theories  are  listed  below  for  reference  and  hereafter  referred 
to  as  Methods  1,  2,  and  3. 


-ms 


-me 


-ms 


Cma 

T 

/s 

00 

Method 

1 

(36) 

cma\ 
.  0  / 

6  sine 

b 

Method 

2 

(37) 

CmaV 
.  0  / 

j  sin20 
b 

Method 

3 

(38) 

Equation  (37)  represents  the  pitching  moment  coefficient  for 
large  angles  (40°)  and  is  an  extension  of  the  linear  form  with 
tail  blanking.  Equation  (38)  is  included  because  it  accounts 
for  the  sin20  variation  of  the  normal  force  coefficient  Cn  which 
has  been  substantiated  for  cylindrical  shapes  for  large  angles. 
The  pitching  moment  slope  coefficient  follows  directly  by  dif¬ 
ferentiating  Equations  (36  ), (37)  and  (38)  giving: 


^mS/00  =  CmCt 


(39) 


9Cms/30  =  !9<Wbsin0 


3Cms/  =  i_(jM  sin2 0  + 

9  0  \  §  /b 


+  (Cma)bcos0 

sin0cos0 
b 


(40) 

(41) 


Equations  (37),  (38),  and  (41)  contain  the  expression  (Cma/0)b, 
which  is  the  nonlinear  variation  of  Cm  with  0  discussed  earl¬ 
ier  . 
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Cms  was  determined  experimentally  for  five  models.  Refer¬ 
ring  to  the  experimental  section  of  the  report  (Appendix  II), 
the  Models  were  3+1,  3+2,  3+3,  3+5,  and  3+6.  The  cylinder  had 
lc/dc  =  2.36  and  the  volute  lv/dv  ranged  from  0.78  to  4.79. 

Cms  for  each  model  along  with  the  theoretical  values  obtained 
from  Equations  (36), (37)  and  (38)  were  plotted  for  each  model 
(Figures  7  through  11).  3+1  shows  the  widest  variation  from 
the  predicted  values  of  the  three  methods.  This  is  so  because 
the  greatest  portion  of  volute  1  (lv/dv  =  0.79)  is  forced  to 
act  in  the  turbulent  slipstream  curling  around  the  cylinder  edge. 
For  this  reason,  short  volutes  should  be  avoided.  Cm(X  varies 
from  a  positive  to  negative  number  throughout  the  first  portion 
of  the  plot  straightening  out  at  around  30  to  40  degrees.  The 
erratic  behavior  is  undesirable  and  should  be  expected  if  lv/dv 
is  less  than  about  1.5  for  dy  -  dc .  The  tendency  toward  a  more 
stable  and  predictable  Cms  with  increasing  lv/dv  is  evident  from 
the  plots.  The  linear  theory  is  surprisingly  accurate  and  can 
be  used  to  predict  Cms .  Method  2  is  more  accurate  and  tends  to 
follow  the  contour  of  the  data  points  for  the  larger  volutes 
(Figure  11) .  Method  3  tends  to  follow  the  data  points  for  0 
larger  than  30  degrees  up  to  the  maximum  range  plotted.  For 
initial  design  it  is  recommended  that  the  linear  theory  be  used. 
If  the  volute  lv/dv  is  less  than  one,  only  wind  tunnel  testing 
should  be  trusted. 

Equations  (40)  and  (41)  are  plotted  (Figures  12  and  13). 

The  Cms  data  was  not  smoothed  and  numerical  differentiation  was 
unproductive  because  the  somewhat  erratic  behavior  of  most  data 
points  was  greatly  exaggerated.  However,  the  Cms  plots  for  the 
five  models  show  that  the  general  slope  of  the  data  points  comes 
very  close  to  the  slopes  of  Methods  1  and  2.  The  closest  agree¬ 
ment  corresponds  to  the  middle  volute  (3) .  The  tabulated  values 
for  Cma/radian  are  as  follows: 


TABLE  I.  Cma/Rad 
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Figure  7.  6  Versus  Cm<;  Model  3+1 
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Figure  12.  6  Versus  Cma  Method  2 
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Since  the  linear  form  is  accurate  enough  for  at  least 
preliminary  design  work,  it  is  possible  to  generate  design 
curves  for  Cma  that  are  simple  to  use.  The  design  curves  pro¬ 
vided  depend  upon  Sv/St,  cnto»  ht >  hbsbcnb0  (designated  B) . 

Recalling  the  definition  of  Cnt0»  it  is  possible  to  plot  Cnt0 
as  a  function  of  Sy/St.  This  is  the  first  design  curve.  Once 
the  forebody  cylinder  has  been  selected  and  volute  stabilizer 
chosen,  Cnt0  can  be  found.  (ht/lo)  is  tabulated  as  a  function 

of  Sv/St  and  Cnt0  in  Table  II.  Since  Sv/St  and  Cnt0  are  now 

known,  (ht/lo)  is  fixed  and  found  by  interpolating  between  the 
two  variables  Sv/St  and  Cnt0*  lo  is  known  since  the  eg  loca¬ 
tion  was  fixed  for  the  selected  f orebody-volute  combination. 

The  tail  parameters,  Cnt0htSt  are  all  known  and  the  forebody 

parameters,  B  =  CnbQhbSb,  are  fixed  once  the  eg  location  is 

specified.  Equation  (8)  can  now  be  solved  for  Cma,  remembering 
that  0  =  0.335.  Equation  (8)  can  be  used  directly  or  the 
master  design  curve  for  Cma  may  be  used. 


(1)  Example:  Find  Gma  for  Model  3+6,  using  Equation  (8) 

Selection  of  volute  6  and  specifying  the  cylinder 
pitch  axis  fixes  the  following  variables: 

Sb  -  2 . 357in2 ;  hb  =  0.81  in;  Cnbo  =  0.833; 

10  =  1.820  in;  St  =  9.027  in2;  Sv  =  6.370  in 
For  Sv/St~  0.7056,  Figure  14  gives  CntQ  =  1.020 

The  ht/l0  matrix  shows: 


Sy / Sf  ' 


Cnto  1.008  1.033 


0.700 

3.264 

3.185 

0.800 

5.619 

5.483 

Interpolating  between  the  values  gives  ht/l0  =  3.226 
Thus,  ht  =  (1 .820) (3.226)  =  5.871  in 
The  parameter  B  =  hbGnbo^b  =  (0.81)  (0 . 833)  (2 . 357)  = 
1.590  in3.  The  quantity  dS  in  Equation  (8)  is 
(1.47) (11,392)  =  16.746  in3 

Inserting  the  numbers  into  Equation  (8) : 

Cma  =-<L5.?.5  ((1.02)  (5.871)  (9 .027)  -  1.59)  =-l. 049/rad 
16.746 
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Figure  14.  Sv/St  Versus  Cnt0 
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Comparing  this  answer  with  Method  1  in  Table  I  shows  a  differ¬ 
ence  of  about  8  percent.  This  follows  from  the  earlier  comments 
regarding  approximating  the  volute  planform  area  as  a  triangle, 
which  was  done  in  the  above  calculation.  Table  II  requires  the 
planform  to  be  triangular  since  this  approximation  was  used  to 
get  ht/l0  in  the  nondimensional  form.  Method  1  answers  in 
Table  I  were  obtained  from  Equation  (33) ,  using  the  more 
correct  trapezoidal  planform  area.  Note  for  volute  6,  dv/dy  = 
4.56. 


(2)  Example:  Find  Cma  for  Model  3+3,  using  the  design 
chart  (Figure  15) . 

The  following  variables  are  fixed  once  again: 

Sb  =  2.357  in2;  hb  *  0.810  in;  Cnb0  =  0.833; 
lo  =  1.820  in;  St  =  4.946  in2  ;  Sv  =  2.289  in2, 

Sv/St  =  0.463 

For  Sv/St  =  0.460,  Figure  14  gives  Cnto  =  0.955 
Find  ht  using  Table  II:  ht  =  2.523  in. 

The  quantity  B  is  the  same  as  example  1,  B  =  1.590  in3 
The  quantity  dS  is:  1.46(7.679)  =  11.21  in3 
Find  the  vertical  axis  labeled  St  and  enter  there 
with  the  value  for  St  =  4.946.  Project  a  line 
and  find  the  line  labeled  ht  =  2.520.  Project  a  line 
perpendicular  to  the  first  through  the  point  defined 
by  the  first  projection  and  the  line  ht  -  2.520.  In 
a  similar  manner,  locate  the  line  labeled  0.958  and 
project  a  third  line  perpendicular  to  the  second, 
and  through  the  point  defined  by  Cnt0  =  0.958  and 
the  second  projection.  Proceed  around  the  graph, 
choosing  next  the  proper  line  for  B  =  1.59  and,  • 
finally,  the  value  for  dS  =  11.21.  Read  the  answer 
on  the  left  side  of  the  vertical  axis  labeled  Cma/9 

Cma/0  =  0.900 

Cma  s  0.302/rad 

Comparison  with  the  answer  in  Table  I  indicates  that  the  pre¬ 
ceding  method  is  high.  This  is  due  to  the  volute  planform  area 
used  in  developing  the  graph  being  triangular.  One  advantage 
to  using  Figure  15  instead  of  Equation  (8)  or  Equation  (33)  is 
that  it  enables  the  designer  to  pick  different  values  for  the 
many  variables  and  see  directly  what  the  effect  upon  Cmg  will 
be.  For  accurate  determination  of  Cm(X,  use  Equation  (33). 

(3)  Example:  Use  Equation  (3  3)  to  find  Cm(X  for  model 
illustrated  below  and  investigate  its  stability. 
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Figure  15.  Master  Design  Curve 


460 


Sb=  2.920  in2;  Sb  -  2.0981  in2;  Sv  =  1.568  in2; 
hb=  1.000  in;  h£  =  0.719  in;  =  2.652  in 
which  gives : 

C  .  .(3, 3811).(0 ,3351.  -  0.118/rad 
“  9.62 

The  low  value  indicates  possible  static  stability 
problems.  Invoking  the  stability  criteria  of 
Equation  (23)  : 

htSt  >  0.883/Cnto^bS^ 

6.64  >  2.560 

According  to  Equation  (23)  the  forebody-volute  combi¬ 
nation  is  stable;  however,  the  volute  is  masked  by 
the  forebody  for  the  early  portion  of  the  oscillatory 
motion.  If  the  volute  lv/dv  is  less  than  1.5,  the 
nonlinear  stability  criteria  stated  by  Equation  (34) 
should  be  used.  In  the  above  case,  lv/dv  =  1.313. 
Equation  (35)  for  the  above  model  gives: 

Cnb^t)^  >  Cn^hbSb,  ®  9  -  0 


and  Cnb0^bSb  =  1*508  <  2.920  =  C^^^hbSb 


This  is  equivalent  to  stating  that  the  model  will  have 
a  trim  angle  different  from  zero.  The  model  will  tend 
to  pitch  over  until  the  volute  balances  the  adverse 
moment  set  up  by  the  large  area  in  front  of  the  eg. 
That  the  model  will  stabilize  out  at  all  was  verified 
by  satisfying  Equation  (23).  The  trim  angle  can  be 
found  by  using  the  follpwing  procedure: 
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0 


CnVQktSv  "  cnb0^bsb  "  Cnb0kbsb>  i.e.,  Cm 


therefore , 


Sv  2  S;  =  0.4046 


Sv  =  1.268  -  0.1839 


(tan6  +  0 . 3808) 2 

tar.9  =  0.0806 
9  =  4°  36 

To  find  the  trim  coefficient  note, 

Cm  =  cmtrim  +  Cma9  =  0 
cmtrim  =  "cma8  =  0.118(0.0785)  =  0.0093/rad 


The  model  was  placed  in  the  wind  tunnel  operating  at  a  dynamic 
pressure  of  about  60  lbs/ft2.  It  immediately  diverged  and 
oscillated  weakly  about  the  10  degree  radial  with  an  amplitude 
ranging  over  approximately  8  to  15  degrees.  A  trim  angle  will 
always  be  present  if  Equation  (35)  is  not  satisfied.  If  this 
is  the  case,  the  magnitude  of  the  angle  will  depend  upon  the 
length  of  the  volute,  small  angles  corresponding  to  long  vol¬ 
utes,  and  large  angles  corresponding  to  short  volutes.  Since 
it  may  be  necessary  ...that  the  forebody  being  stabilized  not 
exceeds  a  certain  trim  angle  after  transients  have  damped  out, 
it  is ' important  to  check  any  design  for  steady  state  trim 
angle.  In  this  respect,  the  above  model  was  analyzed  for 
different  trim  angles  by  varying  the  volute  length.  The  results 
are  shown  in  Figure’ 16. 

Another  point  which  demands  consideration  deals  with  the 
ratio  dv/dc  assumed  to  be  close  to  one.  Actual  fabrication  of 
the  forebody  does  not  allow  this  since  the  volute  must  fit  in¬ 
side  the  cylinder  rim  (Figure  1)  and  thus  dv  <  dc.  This  is 
undesirable  for  stability  since  the  volute  must  now  be  rotated 
a  small  angle  before  being  impinged  upon  by  the  free  stream. 

This  is  distinct  from  the  previous  cases  which  assume  that  the 
volute  starts  producing  a  stabilizing  moment  for  the  slightest 
pitch  of  the  forebody.  The  diameter  of  the  cylinder  and  volute 
must  be  kept  as  close  as  possible,  consistent  with  structural 
requirements.  Test  work  on  a  special  rigid  model  shows  that 
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limit  cycle  vibrations  persist  for  dv  <  0.8dc.  For  larger 
variations,  wind  tunnel  testing  should  be  conducted  to  deter¬ 
mine  if  the  limit  cycles  are  severe  enough  to  warrant  rede¬ 
sign  of  the  volute.  If  the  diamete:.,  cannot  be  brought  close 
enough,  a  longer  volute  can  be  used.  This  has  the  effect  of 
forcing  the  volute  into  the  windstream  early,  reducing  the 
amplitude  of  the  limit  cycle. 

The  damping  derivative  can  be  obtained  from  Equation  (9) 
or  Equation  (42) 


+  cmoP  = 


Cnv0ht,svxvt 


d2S 


1  +  Cnb0sbhbxb^Cnb0Sbh<IjXb 
-  cnVosvhtxvtCnvoSvh^vJ 


(42) 


Using  Equation  (9)  is  less  accurate  because  of  the  triangular 
area  assumption  but  is  quicker  to  use  than  Equation  (42). 
Equation  (9)  yields  the  following  results: 


TABLE  III.  (Cmq  +  Cm£) /Rad (THEORETICAL) 


1+1 

-0.298 

2  +  1 

-0.348 

3  +  1 

-0.506 

1+2 

-0.590 

2+2 

-0.654 

3  +  2 

-0.844 

1  +  3 

-1.124 

2  +  3 

-1.210 

3  +  3 

-1.434 

1  +  4 

-1.666 

2  +  4 

-1.688 

3  +  4 

-1.950 

1  +  5 

-2.476 

2  +  5 

-2.528 

3  +  5 

-2.862 

1  +  6 

-8.048 

2  +  6 

-8.234 

3  +  6 

-8.714 

These  are  plotted  in  Figure  17. 


As  indicated  by  Table  III,  the  damping  coefficient  increases 
dramatically  with  increase  in  tail  length.  Some  experimental 
data  was  obtained  for  Models  3+1,  3+2,  3+3,  3+5  and  3+6. 
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Figure  17.  lv  Versus  (Cmq  +  C,^) 
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TABLE  IV  (Cmq  +  Cmd)/Rad. 


Model 

lv/ dy 

(Cmq  +  Cm^) 

(Cmq  +  Cjtiq)  data 

3+1 

0.785 

-0.506 

-2.370 

0.313 

-0.844 

-2.518 

3  +  3 

1.871 

-1.434 

-3.410 

3  +  5 

2.637 

-2.863 

-4.876 

3  +  6 

4.795 

-8.714 

-10.215 

The  theory  consistently  predicts  low  values  for  (Cma  +  Cm^) .  As 
the  volute  becomes  longer,  agreement  is  better.  Short  volutes 
are  the  greatest  problem  because  forebody  turbulence  is  a  prom¬ 
inent  factor  in  determining  the  flow  oattern  about  the  majority 
of  the  volute  length.  Longer  volutes  operate  at  greater  dis¬ 
tances  and  consequently  are  able  to  act  in  a  more  uniform  flow 
field  of  the  type  the  theory  is  based  upon.  A  second  reason 
for  tlje  large  difference  in  predicted  versus  actual  data  deals 
with  .  In  the  calculation  of  Cma,  ht  had  to  be  averaged  and 
represented  the  only  parameter  that  had  to  be  approximated 
with  average  values. (Cma  +  Cm^)  depends  upon  ht  and  xt  so  that 
averaging  errors  in  both  are  amplified  when  they  are  multiplied 
(Cmq  +  Cm^)  will  have  about  twice  the  error  as  that  in  Cma- 
(Cmq  +  Cm^)  were  obtained  only  for  Models  3+1,  3+2,  3+3,  3+5, 
and43+6.  These  data  points  are  plotted  along  with  the  theoret¬ 
ical  values  for  cylinder  3  models  in  Figure  18. 
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SECTION  III 


FLEXIBLE  BODY  ANALYSIS 


The  various  important  geometrical  parameters  which  affect 
dynamic  and  static  stability  discussed  in  Section  II  provided 
answers  sufficient  for  preliminary  design  work  and,  in  some 
cases  (long  volutes),  good  enough  for  final  design. 

Flexible  volutes  are  not  as  efficient  in  stabilizing 
forebodies  since  they  tend  to  bend  out  of  the  windstream.  In 
addition,  some  configurations  may  be  dynamically  unstable  de¬ 
pending  upon  the  eg  and  the  center  of  pressure  location. 

Static  stability  is  easily  examined  with  slight  modification 
to  the  rigid  body  analysis.  Dynamic  stability  requires  examin¬ 
ing  the  characteristic  equation  of  the  model.  If  the  model  is 
idealized  and  allowed  to  have  only  two  degrees  of  angular  free¬ 
dom,  a  fourth  degree  equation  results  which  can  be  investigated 
with  che  Routh  Hurwitz  criteria.  One  simple  approximation  can 
be  made  to  the  two  equations  of  motion  which  effectively  un¬ 
couple  the  system  and  allow  a  single  equation  to  be  solved. 

This  case  will  be  dealt  with  presently. 

Most  volutes  will  be  light  in  relation  to  the  forebody. 

This  fact  can  be  put  to  use  in  simplifying  the  complex  equations 
of  motion.  If  it  is  also  assumed  that  damping  provided  by 
volute  motion  relative  to  the  forebody  is  small  compared  to  the 
normal  force  set  up  by  the  freestream  velocity,  the  following 
relation  must  hold: 


RCnv  0 . 661ySv$ 

6  =  -  — , - 2 -  8  =  -je 

kv  +  0 . 6  6 1  ySy $ 


(43) 


The  method  using  Equation  (43)  will  be  referred  to  as  the  "J" 
method.  With  Equation  (43),  6  can  be  replaced  and  the  equation 
of  motion  reduced  to  one  variable,  0,  and  the  associated  deriva¬ 
tives.  The  solution  is  not  quite  straightforward,  however, 
because  of  the  proportionality  constant  $.  $  is  found  by  in¬ 

serting  Equation(43)  into  the  $  equation.  Thii,  gives: 


f°r >  $  =  (9+5)  average 

$*0(1  -J)  (44) 
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J,  however,  is  a  function  of  $  so  that  the  above  can  be  written: 


V<P2  +  kv4  -  0kv=  0 


(45) 


where 


V  -  qCnVc)0  •  661ySv 


Note  that  V  can  be  a  function  of  0  depending  on  whether  tail 
interference  is  taken  into  consideration.  If  it  is  not, 

Equation  (45)  can  be  solved  directly  for  a  given  V  and  ky. 

With  this  value,  J  is  known  and  the  problem  is  nearly  solved. 
Since  there  are  no  inertia  effects,  the  aerodynamic  coefficients 
are  changed  by  a  constant  value,  the  new  form  for  each  being 
approximately : 


+  C 


md 


/rigid 


id) 


Cmf /rigid  J 


(47) 


The  two  forms  above  indicate  that  for  flexible  volutes,  the 
overall  static  and  dynamic  stability  decrease  with  increases  in 
Cm<s  and  Cm ; .  This  assumes  that  the  lateral  spring  constant  re¬ 
mains  constant.  Note  that  for  a  given  spring  constant,  in¬ 
creases  in  the  volute  size  cause  .an  increase  in  J  along  with  Cmct, 
cma>  ^ma  changing  the  overall  stability  quickly.  The  conditions 
for  dynamic  and  static  stability  are: 


AB  +VsVCnVod  -  J)  >  0  (48) 

AB*  +  ht'xVtSvCnv  (1  -  J)  >  0  (49) 

L  Vq 

The  stabilizing  effectiveness  of  many  cylinders  is  small  so  that 
in  some  cases  stability  can  be  inferred  by  requiring  that  J  <  1. 
If  the  cylinder  has  positive  AB  and  AB)  then  the  combination 
will  never  be  dynamically  or  statically  unstable,  since  J  <  1 
for  kv  >  0 . 

Thus  far,  the  flexible  analysis  has  not  considered  fore¬ 
body  blanking.  The  first  flexible  forms  are  analogous  to  the 
first  forms  derived  under  the  linear  rigid  body  analysis.  As 
was  shown,  forebody  interference  is  important  and  gives  rise 
to  a  trim  angle.  Returning  to  Equation  (45), V  is  a  function 
of  0  because  Sv  is.  Sv  was  derived  earlier  and  plotted  as  a 
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function  of  6.  To  use  the  results  of  this  work  for  flexible 
models,  it  is  necessary  to  enter  the  Sv  curves  with  the  angle 
(9  -  6).  The  effective  area  is  a  function  of  the  angle  of 
attack  measured  relative  to  the  volute  axis  which  moves  rela¬ 
tive  to  the  forebody  for* flexible  volutes.  To  find  the  result¬ 
ing  pitching  moment  and  its  slope,  it  is  first  necessary  to 
find  J  as  before.  For  a  given  0,  the  deflection  angle  6  of 
the  volute  is  unknown  since  this  is  the  problem.  However,  it 
is  not  possible  to  specify  Sv  without  knowing  <5.  The  calcula¬ 
tion  must  be  an  iterative  process.  Choose  8  and  assume  6  to 
be  zero.  Solve  Equation  (45)  for  9  to  give  the  second  approx¬ 
imation  for  Sv  and  the  procedure  repeats  until  6  approaches  its 
limit  value  and  a  second  angle  chosen,  An  alternate  method  is 
to  assume  what  the  angle  of  attack  (0  -  6)  is,  solve  the  equa¬ 
tions  directly  and  compute  what  cylinder  pitch  angle  will  cause 
the  specified  deflection.  Several  examples  will  illustrate  the 
linear  and  nonlinear  flexible  calculations. 

(4)  Example:  Find  the  pitching  moment  slope  coefficient 
and  damping  coefficient  for  Model  (3  +  6)',  using  the 
linear  method.  Model  (3  +  6)  was  characterized  in 
Example  1.  The  rigid  Cma  is  from  Table  I: 


Cma  =  -0.926/rad 

A  small  leaf  spring  was  inserted  between  the  cylinder 
and  volute,  thus  converting  the  rigid  Model  (3  +  6) 
to  a  flexible  Model  (3  +  6)'.  The  lateral  spring 
constant  was  found  to  be  kv  =  0.27  lbs/rad.  With 
this  new  value,  V  can  be  found: 

V  =  1.124  ft-lbs 


This  form  of  V  assumes  no  interference  effects.  With  V,  $  can 
be  found,  using  Equation  (45). 


$  =  0.188 


J  can  now  be  found 


J  =  0.439 

Cmfi  and  Cm£  are  defined  to  be  (Appendix  I) : 
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-  Cm<$ 


CnynSyht:7 

dcs 


4 


(50) 


-C 


m«5 


cnVos^^t^v 


dcS 


(51) 


Which  compute: 


Cm5  =  -0.510/rad. 
Cm(j  =  -1.614/rad. 


The  aerodynamic  coefficients  follow  from  Equations  (46)and  (47) 


Cmaf=  -0.296/rad. 

:mq  +  cm^£  =  -4.183/rad, 


(5)  Example:  Suppose  that  the  volute  changes  in  length 
in  Example  4,  what  is  the  effect  on  Cm  _?  Table  V 
lists  the  problem  variables  for  1  =  0,  1.118,  1.825, 

2.62,  3.89,  7  and  «>. 


TABLE  V.  Cnw/Rad 


lv 

V 

4 

j 

cma 

Cm5 

cmaf 

0 

0 

A 

0 

0 

-0.0558 

0 

-0.0558 

1.118 

0.0281 

0.324 

0.0326 

-0.126 

-0.103 

-0.119 

1.825 

0.0721 

0.309 

0.0763 

-0.196 

-0.167 

-0.168 

2.62 

0.288 

0.1403 

-0.275 

-0.230 

-0.2024 

3.89 

0.343 

0.245 

-0.470 

-0.344 

-0.2707 

7.00 

1.124 

0.439 

-0.926 

-0.510 

-0.2958 

CO 

OO 

0 

1 

-  CO 

0 

'vO 
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As  indicated  by  Table  V,  there  is  a  limit  to  the  effectiveness 
of  a  flexible  volute  beyond  which  little  is  gained  with  fur¬ 
ther  increase  because  of  the  tendency  to  bend  out  of  the  wind- 
stream.  This  is  distinct  from  the  rigid  case  where  increase 
in  volute  length  is  manifest  in  greater  Cma  without  limit. 

Larger  flexible  volutes  provide  greater  righting  moments  than 
smaller  volutes;  however,  their  efficiency  decreases  as  indi¬ 
cated  by  a  falling  off  of  Cma^  beyond  a  certain  best  length  ly. 
If  nothing  else  mattered,  Equation  (46)  would  be  maximized  and 
that  value  of  lv  used.  As  lv  =  5.5  inches,  Cma  reaches  about 
-0.290  and  an  increase  of  lv  by  1.5  to  7  inches  increases  Cmar 
by  1.8  percent.  In  the  limit,  lv  *  °°>  J  goes  to  1,  $  to  0  ana 

Cm  to  approximately  0.  To  find  the  best  design  length,  the 
af 

following  example  furnishes  complete  calculations. 

(6)  Example:  Show  the  effect  of  increasing  ly  for  a 
cylinder  with  lc/dc  =  2  where  xCg  =  1.46  inches  and 
dc  =  1.46  inches.  Table  VI  shows  the  problem  values. 

A  plot  of  Table  VI  gives  Cmaf  a  maximum  value  of  Cmaf  = 
-0.3007  at  lv  =  7.5  inches.  Beyond  this,  a  larger 
volute  would  reduce  in  effectiveness.  The  curve 
approaches  its  maximum  value  rapidly  and  this  value 
is  nearly  attained  with  smaller  volutes.  If  the 
best  design  value  is  defined  to  be  90  percent  (x)lymax 
(maximum) ,  a  more  reasonable  length  volute  can  be 
used.  In  this  case,  the  best  design  value  is  -0.2764, 
which  occurs  at  about  lv  =  4.5  inches  and  represents 
a  77  percent  reduction  over  the  maximum  value  volute 
length  (Figure  19) . 

All  of  the  foregoing  deal  with  the  linear  form  of  the 
flexible  equations;  no  account  was  made  for  forebody  inter¬ 
ference.  Thus,  the  area  Sv  used  for  all  calculations  was  the 
entire  volute  planform  area.  To  account  for  interference  and 
find  the  value  of  0tfim>  it  is  necessary  to  invoke  the  results 
of  Equation  (29).  As  suggested  earlier,  the  procedure  is  itera¬ 
tive  . 

(7)  Example:  Find  Cmct  for  Flexible  Model  1,  using 
iteration.  The  calculation  will  he  carried  out  for 
0  =  0°,  5°,  10°,  15°,  20°,  25°,  30°,  35°  and  40°. 
Assume  6=0  and  9=5°. 

For  the  ratio  lv/dc  =  2.915,  Equation  (29)  gives 
Sy/a-2/  =  0.868  and  Sy  =  1.1487  in.2  With  this 
number,  the  problem  can  be  solved  as  in  Example  4. 

This  gives  J  =  0.1012.  Thus,  5  -was  not  zero  but  is: 

6  =  -  J9  =  -  (0.1012)5°  =  0.506° 
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TABLE  VI.  BEST  DESIGN  LENGTH  COMPUTATIONS 


O 

LO 

rH 

to 

H- 

rH 

to 

CM 

H* 

© 

00 

CTl 

LO 

© 

m 

H- 

CM 

L0 

to 

\o 

LO 

On 

o 

00 

to 

© 

o 

o 

00 

3 

to 

rH 

LO 

oo 

rH 

to 

LO 

nO 

oo 

00 

cr> 

© 

© 

o 

o 

© 

s 

o 

rH 

rH 

rH 

(M 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

to 

to 

CM 

U 

• 

• 

• 

• 

. 

• 

* 

« 

• 

• 

• 

• 

t 

• 

• 

• 

o 

O 

1 

O 

l 

o 

1 

O 

i 

o 

1 

o 

O 

l 

o 

o 

1 

o 

1 

O 

1 

O 

1 

o 

1 

o 

1 

o 

1 

O 

vO 

00 

on 

LO 

00 

LO 

o 

rH 

o 

nO 

r-* 

00 

CM 

CM 

© 

LO 

© 

CM 

co 

LO 

rH 

nO 

On 

© 

00 

LO 

rH 

LO 

[n 

00 

o 

'S’ 

rH 

oo 

CM 

LO 

oo 

rH 

H* 

O 

CM 

H- 

© 

o 

o 

o 

rH 

rH 

rH 

CM 

CM 

CM 

to 

to 

to 

H* 

M- 

o 

o 

o 

o 

O 

° 

O 

O 

o 

O 

o 

o 

O 

o 

O 

o 

to 

On 

to 

rH 

00 

CO 

CM 

CM 

CT> 

iH 

rH 

© 

CM 

LO 

to 

cn 

cO 

r- 

LO 

CM 

o 

CTi 

00 

00 

© 

o 

CM 

LO 

00 

to 

to 

rH 

o 

On 

CO 

nO 

rt 

to 

CM 

rH 

O 

o 

© 

O0 

c^. 

to 

to 

to 

to 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

rH 

rH 

o 

o 

o 

o 

O 

o 

O 

O 

O 

© 

O 

O 

O 

o 

o 

O 

O 

o 

\0 

rf 

© 

o 

vO 

H- 

nO 

o 

vO 

LO 

L0 

H- 

rH 

to 

LO 

LO 

LO 

CM 

o 

nO 

CM 

r-> 

rH 

•<* 

00 

© 

> 

o 

H’ 

rH 

nO 

CM 

On 

no 

LO 

© 

r- 

o 

to 

CM 

o 

o 

° 

rH 

rH 

CM 

CM 

to 

H- 

LO 

© 

h- 

© 

o 

rH 

to 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

rH 

rH 

rH 

o 

00 

o 

GO 

LO 

LO 

CTi 

o 

CT> 

to 

CM 

© 

a 

vO 

o 

to 

H- 

o 

LO 

CO 

O 

o 

rH 

to 

rH 

o 

to 

to 

CM 

<2 

o 

CM 

© 

o 

LO 

o\ 

CM 

nO 

© 

CM 

LO 

t-- 

O 

CM 

LO 

u 

to 

rH 

rH 

CM 

CM 

CM 

to 

to 

to 

H- 

L0 

LO 

LO 

H* 

o 

1 

O 

O 

l 

o 

I 

O 

1 

o 

1 

o 

i 

o 

o 

O 

O 

t 

o 

i 

o 

1 

o 

t 

O 

1 

O 

t 

o 

o 

to 

on 

to 

nO 

00 

vO 

CM 

rH 

CM 

LO 

© 

to 

to 

t". 

3 

LO 

© 

00 

LO 

nO 

o 

nO 

rf 

to 

to 

H* 

© 

© 

to 

t-- 

s 

to 

CM 

to 

On 

LO 

CM 

00 

LO 

© 

© 

to 

o 

00 

LO 

u 

o 

rH 

rH 

CM 

CM 

to 

H" 

LO 

no 

© 

00 

* 

© 

o 

o 

1 

O 

l 

o 

1 

O 

i 

O 

i 

o 

1 

O 

* 

o 

1 

O 

1 

o 

1 

O 

l 

o 

i 

o 

1 

o 

1 

o 

1 

o 

1 

H- 

to 

© 

00 

CM 

t^- 

o 

to 

nO 

© 

CM 

LO 

o 

rH 

to 

cm 

LO 

CM 

LO 

00 

CM 

h* 

00 

CM 

LO 

00 

iH 

LO 

00 

rH 

LO 

00 

o 

OJ 

t"- 

00 

to 

00 

On 

•H- 

o 

LO 

o 

© 

rH 

© 

CM 

r^. 

CM 

TJ 

• 

• 

• 

• 

♦ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

© 

© 

r-~ 

00 

00 

on 

On 

o 

rH 

rH 

CM 

CM 

to 

to 

M- 

H* 

LO 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

o 

LO 

LO 

o 

LO 

LO 

LO 

o 

L0 

o 

LO 

O 

LO 

> 

© 

a\ 

VO 

CM 

O) 

LO 

CM 

00 

LO 

rH 

CO 

H* 

rH 

o 

CO 

to 

o 

00 

rH 

LO 

on 

CM 

VO 

O 

to 

rH 

■Sf 

00 

CM 

o 

rH 

O 

rH 

CM 

CM 

CM 

to 

to 

H* 

LO 

LO 

LO 

© 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

\£> 

VO 

vO 

O 

nO 

vO 

nO 

nO 

nO 

vO 

no 

© 

vO 

© 

© 

© 

© 

pa 

on 

o\ 

a> 

CT> 

cr> 

On 

Oi 

on 

© 

© 

© 

© 

© 

© 

© 

© 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

vO 

’O 

vO 

vO 

\o 

nO 

vO 

vO 

nO 

vO 

© 

© 

© 

v£> 

vO 

© 

© 

on 

on 

On 

on 

on 

CT> 

o> 

cn 

CT> 

CT> 

© 

© 

a* 

© 

© 

© 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

-  -P 

© 

vO 

LO 

00 

rH 

H- 

O 

to 

VO 

© 

CM 

LO 

00 

rH 

H- 

r-« 

<3* 

rH 

H* 

rH 

H* 

o 

H* 

O 

* 

t"- 

o 

rH 

rH 

CM 

CM 

to 

to 

to 

H- 

L0 

LO 

LO 

vO 

vO 

© 

> 

O 

LO 

LO 

O 

LO 

o 

LO 

o 

LO 

O 

LO 

o 

L0 

O 

LO 

O 

LO 

rH 

♦ 

* 

• 

• 

• 

• 

• 

• 

• 

• 

* 

• 

• 

• 

• 

o 

rH 

CM 

CM 

to 

to 

LO 

LO 

© 

vO 

Is- 

00 

00 

42 


f 


I 


I 

V 


8 


t/i 

3 

in 


>  CJ1 


<u 

)-< 

3 

W> 

•H 

fcU 


43 


i 


(Inches) 


The  net  angle  of  attack  is: 

6volute  =  Beyl."  6  =  5°  -0.506°  =  4.494° 

Enter  Equation  f 29 )  with  this  new  angle  and  find 
Sv  =  1.081  in.2  Repeat  the  procedure.  The  following 
table  of  values  results: 


TABLE  VII.  Cmaf/Rad. 


0 

Sv 

<J> 

J 

6 

cm<$ 

cma 

cmaf 

0° 

0 

0.335 

0 

0 

0 

0 

0 

5° 

1.0889 

0.3026 

0.0968 

-0.4413 

-0.1642 

-0.1799 

-0.1466 

O 

o 

H 

1.484 

0.2938 

0.1238 

-1.102 

-0.2172 

-0.2492 

-0.1917 

na 

1.677 

0.2894 

0.1361 

-0.2418 

-0.2856 

-0.2’39 

ISM 

1.786 

0.2871 

0.1430 

-2.502 

-0.2555 

-0.3055 

-0.2236 

25° 

1.85S 

0.2857 

0.1475 

-3.204 

-0.2640 

-0.3195 

30° 

1.903 

0.2848 

0.1498 

-3.890 

-0.2700 

-0.3266 

-0.2372 

35° 

1.936 

0.2844 

0.1517 

-4.592 

-0.2743 

-0.3563 

-0.2439 

40° 

1.959 

0.2837 

0.1531 

-5.308 

-0.3048 

-0.3397 

-0.2410 

The  results  of  the  flexible  test  work  for  flexible  Models  1 
and  2  are  plotted  and  the  predicted  values  are  plotted  (Figures 
20  and  21).  As  was  done  in  the  rigid  analysis,  Cm  was 
chosen  to  represent  the  pitching  moment  data  and  tne  three 
analytical  curves  represent  the  same  three  theories  outlined  in 
the  rigid  work.  The  first  is  the  most  linear  flexible  theory 
('j.e.,  Equation  (46)  -  no  forebody  blanking)  and  the  second 
and  third  take  into  account  $sin  (9  +  6)  and  sin2  (0  +  6) 
variations  of  pitching  moment  plus  forebody  blanking. 

As  previously  stated,  the  flexible  volute  operating  be¬ 
hind  the  cylindrical  forebody  will  have  a  trim  angle  different 
from  zero  if  the  eg  of  the  model  is  located  aft  of  the  eg  of 
the  cylinder.  Returning  to  Example  4,  if  the  volute  is  imag¬ 
ined  to  have  elasticity  with  kv  =  0.27,  the  following  example 
can  be  solved  for  the  trim  angle: 

(8)  Example:  For  this  a  trapezoidal  area  will  be  used. 
The  rigid  body  analysis  provided  two  equations 
which  can  be  used  for  this  purpose.  The  total  area 
into  the  windstream  at  a  given  angle  of  attack  is 
the  sum  of  the  following  two  expressions: 
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Pitch  Angle  6  (Degrees) 

Figure  21.  0  Versus  Cmj_  Flexible  Model 
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Sy"  =  dy  ( ly  -  X  ) 


Thus 


Sv  = 


V 


+ 


s 


ti 

v 


The  condition  for  finding  the  trim  angle  is  found  by  setting 
Equation  (48)  to  zero. 


AB  +  htSvCny0(l  -  J)  =  0 


Recalling  the  values  of  the  constant  terms  from  Example  4,  the 
above  can  be  written  as: 


1.568  -  0.186  -  0.0818 _  =0.4046 

tan9  +0.291  (tan0  +0.  2907)  2  l  -  J 


An  iterative  procedure  may  be  used.  The  following  table  gives 
intermediate  answers  and  the  final  value  for  the  trim  angle. 


TABLE  VIII.  TRIM  ANGLE 


9trim 

Sv 

V 

$ 

J 

Left 

Side 

Right 

Side 

0 

0 

0 

0 

0 

-0.039 

0.4046 

30' 

0.057 

0.005 

0.333 

0.006 

0.037 

0.4070 

1° 

0.109 

0.009 

0.332 

0.009 

0.108 

0.4083 

2° 

0.197 

0.017 

0.328 

0.020 

0.226 

0.4129 

4° 

0.336 

0.029 

_ 

0.324 

_ 

0.034 

0.419 

0.419 

It  remains  to  look  into  the  damping  derivatives,  using 
the  J  method.  These  calculations  are  quite  simple  and  can  be 
carried  out  quickly,  using  Equations  (47)  and  (51) .  This  was 
done  for  the  models  listed  in  Example  5  and  are  tabulated 
below : 
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TABLE  IX.  (Cmq  +  Cm^) £/Rad .  (THEORETICAL) 


Model 

lv 

(Cmq  +  Cm<j) 

$/o 

Cm,$ 

(Cmq  +  Cmd)f 

3  +  1 

1.118 

-0.506 

0.967 

-0.169 

-0.484 

3  +  4 

1.805 

-0.844 

0.922 

-0.315 

-0.754 

3  +  3 

2.62 

-1.434 

0.859 

-0.519 

-1.159 

3  +  5 

3.89 

-2.862 

0.755 

-0.891 

-1.904 

3  +  6 

7.00 

-8.714 

0.561 

-1.871 

-4.067 

The  J  method  represented  a  good  first  approximation  for 
determining  the  effects  of  volute  flexing  upon  the  basic 
stability  derivatives.  The  gi  latest  shortcoming  deals  in  the 
area  of  dynamic  stability;  such  instability  cannot  be  predicted 
since  the  volute  is  assumed  massless.  Wind  tunnel  testing 
shows,  however,  that  even  for  light  volutes  of  the  type  which 
would  probably  be  used,  dynamic  instability  is  a  serious  problem. 

A  coning  motion  which  persists  after  the  initial  tran¬ 
sients  have  damped  out  is  apparently  caused  by  turbulence  and 
coupling  between  the  pitching,  yawing  and  stretching  motions  of 
the  volute.  The  three  motions  seem  to  enhance  each  other,  per¬ 
iodically  damping  out  only  to  be  revived  to  the  same  or  greater 
amplitude  motion.  For  speeds  lower  than  Mach  0.19,  the  motion 
was  less  violent.  The  following  theory  deals  only  with  two  di¬ 
mensional  motion  and  allows  for  only  two-degree  angular  freedom 
motion  in  one  plane. 

Equations  (1)  and  (2)  represent  the  backbone  of  this  theory. 
To  fully  grasp  the  meaning  of  each  leading  coefficient,  the  two 
equations  will  be  broken  down  to  the  individual  parts  and  dis¬ 
cussed  separately.  The  basic  equations  for  reference  are  again: 

a  8  +  a  6  =  -  a  0  -  a  9  -  a  <5  -  a  6  (52) 

’-1  12  13  14  15  16 

a  9  +  a  6  =  -  a  0  -  a  0  -  a  5  -  a  6  (53) 

21  22  23  24  25  26 

a  -  The  inertia  of  the  forebody  and  volute  are  accounted  for 
11  with  aM.  This  term  is  the  same  as  that  which  would 
appear  with  8  for  a  rigid  model. 


=  Ib+  Iv  +  2  MvRvFy  +  MVRV2-  I  (54) 
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ai2  -  Allowing  the  volute  to  oscillate  and  giving  it  mass  and, 
hence,  inertia,  gives  rise  to  an  inertia  force  which  is 
transmitted  to  the  forebody  through  the  elastic  coupl¬ 
ing.  This  inertia  force  must  be  included  on^the  inertia 
side  of  the  moment  summation  equation.  For  6=0.  The 
a  inertia  term  drops  out  of  the  equatior  and  the  left 
siae  of  Equation  (52)  reduces  to  the  rigid  form. 


al2  ~  iv  +  ‘V  Rv  ty  =  Ky.  (55) 

a13  -  lhis  term  represents  the  damping  moment  set  up  by  a  rigid 
body  model  having  the  same  tail  and  forebody  as  the  flexi 
ble  one.  The  term  is  relatively  large  and  represents  a 
major  contribution  to  the  overall  damping  of  the  model. 
The  term  has  an  order  of  magnitude  of  about  10' 3 .  It  can 
be  written  as: 


dc<i> 

-  (56) 


or,  using  linear  aerodynamic  notation: 


Tpv-  Sd( 


^^nb0^bxb^b  2CnkoS]/  xt>  hjj  2CnVoSvxVthi: 

J 


dr  2  S 


AB' 


dr  2S 


dr2S 


-qdcS 


(57) 


(Cm  +  cm • )  is  the  same  as  that  defined  for  the  rigid 
modSls.  An  important  fact  yet  to  be  shown  is  that  a13 
must  be  large  m  relation  to  ai 5  [its  complimentary  term 
in  Equation  (52)  ]to  aid  in  insuring  dynamic  stability. 

The  terms  in  Equation  (56)  indicate  how  to  make  ai  3  large 
This  can  be  accomplished  without  increasing  a15  by  in¬ 
creasing  the  forgbody  size  or  by  moving  the  eg  forward 
which  increases  xyt.  This  latter  point  will  be  examined 
more  thoroughly  since  the  same  adjustment  will  provide 
better  static  stability.  xVt  plays  a  critical  role  in 
determining  whether  the  model  will  have  tendencies  toward 
dynamic  instability. 


a14  -  This  is  simply  the  rigid  expression  for  static  pitching 

moment.  It  is  one  of  the  largest  terms  in  Equation  (52), 
having  an  order  of  magnitude  of  about  10' l.  The  term 
must  be  positive  if  the  model  is  to  be  statically  stable. 
It  can  be  written  two  ways  as: 
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2 


(58) 


a 


1  5 


-  V 

2  p  o 


or 


'nb0^bhb  +  Cnj)oSb/  hb'  +  CnVQSvht 
dpS  dpS  dpS 

—  qdpS  Cihqj 


(59) 


This  term  does  not  show  up  favorably  since  it  causes 
trouble  in  the  characteristic  equation  of  the  model. 

A  sign  reversal  on  one  coefficient  of  the  characteristic 
equation  will  occur  if  a15  is  too  large  which,  according 
to  the  Routh-Hurwitz  criteria,  guarantees  dynamic  in¬ 
stability.  Unfortunately,  a15  cannot  be  indiscriminately 
reduced  since  several  of  the  constituent  terms  are  re¬ 
quired  to  be  large  for  large  a1?  and  a14  which  are  con¬ 
sidered  to  be  favorable  quantities  (i.e.,  the  larger, 
the  better).  a} s  is  written  two  ways  as: 


_i_ 

2P 


V  2Sd 

nr>  ** 


C 


2CnVc)SvXvht/ 

d2pS 


dp$ 


2V0 


(60) 


and,  by  using  linear  aerodynamic  notation,  as: 


-qdcS  (61) 

It  would  be  difficult  to  have  large  at 3  $nd  small  at g 
if  not  for  the  subtle  difference  between  Xy  and  xvt 
because  the  forebody  AB'  is  usually  small.  xvt  can  be 
written  as  a  function  of  the  volute  length  and  the  eg 
location.  For  a  given  volute,  xvt  controlled  by 
the  eg  location.  Movement  forward  will  be  manifest  in 
greater  xvt  which  helps  to  insure  dynamic  stability. 
xv  is  not  affected  by  eg  location  and,  hence,  does  not 
change.  If  the  eg  were  to  shift  backwards  and  longer 
volutes  were  added  to  the  same  forebody,  aj 5  would  in¬ 
crease  at  a  faster  rate  than  a!  3  which  is  undesirable. 
xv  and  xvt  are  both  related  to  center  of  pressure  and 
their  definition  carefully  takes  into  account  the  volute 
planform  shape.  Thus,  the  above  indicates  how  the  center 
of  pressure  and  the  eg  location  need  to  be  related. 

-  This  term  is  analogous  to  a14  and  represents  the  static 
moment  set  up  by  deflection1 of  the  volute  into  the  wind. 
Note  that  this  deflection  is  measured  relative  to  the 
forebody.  Depending  on  whether  the  volute-  bends  into 
the  wind  or  away  from  it,  a  can  add  or  detract  from 
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static  stability.  Like  a  ,  alg  does  not  show  up  favor¬ 
ably  in  the  dynamic  stability  criteria.  If  the  forebody 
has  the  eg  location  in  its  geometrical  center,  then  a16 
is  equal  to  a14.  This  is  evident  from  the  defining 
equation : 


1  V  2 

2  p  «° 


Sdc 


Jn 


Vo 


S„h 


v"t 


L  dcS 


or  as  : 


-qdcS 


Cms 


(62) 

(63) 


Thus  a16  has  the  same  order  of  magnitude  as  a14.  There 
is  no  way  to  increase  a14  without  also  increasing  AB  if 
h£  is  measured  from  the  model  eg  and  is  the  same  for 
both  coefficients. 


The  second  equation  of  motion  was  found  through  a  moment 
summation  about  the  tail  hinge  point.  Because  of  this, 
it  was  necessary  to  alter  the  normal  expression  for  the 
inertia  terms  and  include  the  effects  of  a  linearly 
accelerating  coordinate  system.  When  the  inertia  side 
of  the  moment  summation  is  written  out,  a2l  appears  as 
the  lead  coefficient.  In  terms  of  the  physical  con¬ 
straints  : 


Jv  +  MvRvrv  (64) 

This  is  the  second  portion  of  the  inertia  expression. 

It  is  defined  to  be: 


Iv  (65) 

When  a13  has  been  determined,  a?J  can  be  found  by  sub¬ 
tracting  the  forebody  contribution  to  a.3  and  multi¬ 
plying  by  0.661v/h't.  az3  represents  a  clamping  moment  set 
up  by  forebody  pitching.’  Note  here  that  the  moment  arm 
is  only  0.661v,  since  the  moment  summation  center  was 
located  at  the  hinge  line.  a2 3  can  be  defined  in  several 
alternate  ways,  as  follows: 


2Cyiy^SyXy.£0 . 661y 


dr2S 


2V  _ 


(66) 
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0  »  661y  Xy  dp 

-qdcS  (CH)— 7!  Jt  -1-  (67) 

V  '  ht  xv  2Vro 

/  \0.661v 
[  -  (forebody  contribution)  1  -  (°8) 

a2,  has  a  destabilizing  effect  on  the  dynamics  of  the 
model.  The  same  conclusion  previously  reached  concerning 
eg  versus  center  of  pressure  location  can  be  arrived  at 
here  by  examining  the  moment  arm  Q.661v.  Note  that  a13> 
which  serves  the  same  variable  in  Equation  (52)  as  does 
a23  in  Equation  (53),  depends  upon  h£  which  itself  relies 
on  eg  location  for  a  fixed  volute.  Increasing  rearward 
location  of  the  eg  must  be  done  with  some  care  if  dy¬ 
namic  stability  is  to  be  maintained. 

This  term  is  the  static  pitching  contribution  due  to 
excursions  of  e  in  Equation  (53).  It  has  a  similar  form, 
as  does  a14,and  can  be  written  in  several  ways: 


I  r  ^ 

2p'w 


CnVosv  3li 


-qdcS(Cm6)  0.661v/h£  (70) 

-  (forebody  contribution)^ ‘ ^61y 

• 

Dynamic  stability  is  aided  if  a25  is  large.  One  finds 
here  that  a  compromise  must  be  made  among  the  variables 
since  the  previous  findings  concerning  eg  location 
versus  large  lv  are  not  manifest  in  a25.  That  is,  a25 
can  be  made  large  regardless  of  eg  location.  By  increas¬ 
ing  lv,  however,  to  enlarge  a25,  destabilizing  elements 
are  at  work  through  a15  and  a23,  both  of  which  are  grow¬ 
ing  with  increases  in  lv.  a13  also  increases  with  lv. 
The  relation  between  these  coefficients  to  determine 
whether  such  changes  in  the  model  geometry  cause  insta¬ 
bility  is  the  subject  of  the  Routh  Hurwitz  criteria. 
a25  may  be  written  as: 

1  2Cj}  SyXy0.661y  dr$ 

i0V„2Sdc  - 7  JL  (72) 

L  d2cS  J  2V*, 
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(73) 


1 2  6 


(a 

1  3 


-qdcS 


(cm<$) 


0 . 661v 


xv0 . 661v 

forebody  contribution)  * - — 

xvth  t 


(74) 


-  This  is  a  unique  coefficient  and  includes  the  effects  of 
the  elastic  coupling.  It  is  the  only  coefficient  that 
includes  the  spring  constant  kv.  It  has  already  been 
shown  that  the  value  of  the  spring  constant  is  important 
in  determining  the  value  of  pitching  and  damping  moments 
set  up  by  the  model's  motion.  The  dynamic  criteria  shows 
that  large  kv  helps  to  provide  a  dynamically  stable  model. 
a2,  is  a  function  of  a24  which  can  be  typified  as  an 
unfavorable  value  if  too  large.  For  a26  to  be  large 
without  compromising  stability  through  a24,  ky  must  be 
large.  This  stands  to  reason  since,  in  the  limit  as  l:v 
approaches  infinity,  the  system  reduces  to  its  rigid 
body  form.  a26  is  easiest  to  use  defined  as  below: 


All  of  the  preceding  comments  regarding  dynamic  stability 
result  from  the  Routh  Hurwitz  criteria.  An  absolute  stability 
criteria  involving  the  model  constraints  was  not  obtained  be¬ 
cause  of  the  great  amount  of  algebraic  work  and  the  practically 
impossible  task  of  using  the  results.  Thus,  the  volute  equipped 
body  has  not  been  proven  to  be  unstable  for  a  given  type  of  con¬ 
figuration.  To  prove  that  volute  equipped  models  cannot  be  un¬ 
stable  requires  writing  the  coefficients  of  the  characteristic 
equation  in  the  stability  criteria  scheme  and  insuring  that  the 
proper  terms  will  not  change  sign  by  relating  the  constraints 
back  to  the  original  a^jig. 

It  was  found  by  computing  many  examples  for  a  variety  of 
different  volute  and  forebody  combinations  that  a  good  indi¬ 
cator  of  impending  dynamic  instability  was  the  coefficient  C 
of  the  s3  term  of  the  fourth  degree  characteristic  polynomial. 
This  term  will  be  the  first  one  to  change  sign,  thus  indicating 
dynamic  instability.  That  is  not  to  say  that  no  change  is 
reason  not  to  invoke  the  remainder  of  the  stability  criteria. 

If  the  mentioned  coefficient  is  much  smaller  than  the  other 
coefficients  but  does  not  change  sign  (all  coefficients  would 
then  have  the  same  sign) ,  there  is  a  good  chance  the  model  will 
still  have  two  unstable  roots  and  two  stable  roots.  The  Routh 
Hurwitz  criteria  is  discussed  more  fully  in  Appendix  I.  From 
that  discussion  the  following  observations  are  made  which  tie 
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the  preceding  to  the  a^ig.  The  fourth  degree  characteris¬ 
tic  equation  can  be  written  as: 

C j  j s 4  +  C12s3  +  C13s2  +  C14s  +  C1S  =  0  (76) 

The  coefficient  of  interest  is  C12.  C12  can  be  written  in  terms 
of  the  ajj's  as  follows: 


1  2 


=  (aM. 


'2  5 


a2aai 3^  - 


(a  a 
v  2  1  1 


a  a  ) 

1  2  2  3  ' 


(77) 


Since  CJ2  may  be  the  first  of  the  Cjj’s  to  go  negative,  a 
necessary  but  not  sufficient  condition  for  dynamic  stability 
is : 


aa  +aa>aa  +aa 

1125  2213  2115  1223  (78) 

This  inequality  constraint  is  a  relation  between  the  inertia 
of  the  cylinder  and  volute  and  the  damping  properties  of  the 
cylinder  and  volute  combination.  Note  that  the  spring  con¬ 
stant  does  not  appear  in  Equation  (78).  The  small  mass  assump¬ 
tion  used  in  the  J  method  will  reduce  Equation  (78)  to: 


aa>0,a  >0,  a  >0  (always)  (79) 

1125  11  25 


Thus  massless  volutes  will  probably  be  stable.  Giving  the 
volute  mass  requires  matching  the  ajj's  so  that  C12  cannot 
change  sign.  Note  that  C  will  usually  be  positive  which  re¬ 
quires  that  C  j  2  be  greater  than  zero. 


Equation  (78)  is  useful  in  finding  out  if  a  given  confi¬ 
guration  will  have  an  obvious  dynamic  stability  problem.  Sub¬ 
stituting  for  the  volute  and  cylinder  geometrical  and  inertia 
properties  for  the  a^j's,  the  following  equation  results: 


K 


v 


rigid 

L 


+ 


0.661v 


< 


cmq  + 


(80) 

rigid 


Note  that  Kv,  h^  Iv,  xVt,  lv  are  all  interdependent  and 
changing  any  one  changes  the  others.  Increasing  lv  causes 
the  remaining  parameters  to  increase.  I,  however,  can  be 
changed  without  affecting  the  other  variables  by  increasing  Sb. 
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This  is  a  favorable  trend,  since  Equation  (80)  shows  that  I  re¬ 
duces  the  left  side  and  enhances  the  inequality.  Having  a 
large  cylinder  in  relation  to  the  volute  tends  to  produce 
greater  dynamic  stability,  holding  all  other  parameters  the 
same.  Equation  (80)  can  be  used  to  examine  the  stability  of 
several  cylinder  volute  combinations. 

Table  X  lists  the^aijs  for  several  models.  These  values 
were  calculated  using  xv  =  0.661v  and  xVt  =  10  +  0.5001v. 

Note  that  here  an  approximation  was  made  since  the  earlier 
rigid  work  indicated  that  these  variables  h£  and  hv  will  be 
functions  of  0  and  0.  h-^  was  set  equal  to  10  +  |-lv  and  hv  = 
0.061v.  The  ajis  were  multiplied  by  1000  which  gave  the 
following  set  or  equations: 

Model  (3+1)' 

0.2040*  +  0.036^  +  0.0430  +  380  +  0.00926  +  336  =  0 
0.0360  +  0.0116  +  0.0090  +  9.600  +  0.0027$  +  279.66  =  0 


Model  (3+2) * 


0 . 343  0*  +  0.1126  +  0.0770  +  660  +  0.028$  +  616  =  0 
0.1120  +  0.0586  +  0.0250  +  24.40  +  0.0116  +  294.46  *  0 


Model  (3+3)' 


0.605  0*  +  0.272  6  + 
0.2720  +  0.1566  + 


0.1240  +  1030  +  0.064$  +  986  =  0 
0.0616  +  51.50  +  0.0346  +  321.56  =  0 


Model  (3+5)' 


1.2960  +  0.7375  + 
0.7350  +  0.4606  + 


0.291$  +  2030  +  0.1936  +  1986  =  0 
0.1660  +  116.40  +  0.1156  +  386.46  =  0 


Model  (3+6)’ 


5.7600  +  4.2226  + 
4.2226  +  3.1835  + 


1.0780  +  5350  +  0.9296  +  5306  =  0 
0.7690  +  381.40  +  0.675$  +  651.46  =  0 


From  these  equations,  the  characteristic  equation  can  be  derived 
by  noting  that  the  initial  values  for  0  and  6  will  not  be  zero. 
Substituting  for  the  variables  according  to  the  scheme  in 
Appendix  I,  the  following  set  of  equations  can  be  written: 


Model 

(3  +  1)  ’ 

s4 

+  0.295s3 

+  25,850s2 

+  5,400s  +  4,746,000 

Model 

(3  +  2)  ’ 

s4 

+  0.313s3 

+  12,960s2 

+  2,800s  +  2,440,000 
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Model  (3+3)' 


s 4 

Model 

+  0.289s3 

(3+5)' 

+  8,330s2 

+  1600s  + 

1,376,000  =  0 

s4 

Model 

+  0.335s3 

(3+6)' 

+  6,200s2 

+  1400s  + 

947,000  =  0 

s 4 

+  0.278s3 

+  3,050s2 

+  570s  + 

278,000  =  0 

To  find  the  roots  of  these  equations,  it  is  necessary  to  use  the 
results  of  Appendix  I  which  outline  a  graphical  scheme.  It  is 
helpful  and  also  time  saving  if  the  characteristic  equations  are 
examined  with  the  Routh  Hurwitz  criteria  before  attempting  to 


factor  them. 

The  values  of  the  criteria  terms 

Model  (3+1)' 

1 

0.C  95 
0.011 
0.125 

25,850  4,746,000 

5,400 

6.607 

Model  (3+2)' 

1 

0.313 

0.067 

1.295 

12,960  2,440.000 

2,800 

41.32 

Model  (3+3) » 

1 

0.289 

0.335 

10.026 

8,330  1,376,000 

1600 

165.991 

Model  (3+5)' 

1 

0.335 
7.038  2 

458.451 

6,200  947,000 

1400 

,443.452 

Model  (3+6)’ 

1  3,050  278,000 

0.278  570 

65.074  12,392.047 

18,542.087 
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Since  the  first  column  of  each  scheme  does  not  change  sign,  all 
given  models  are  dynamically  stable.  Thus,  the  quadratic  fac¬ 
tors  of  the  fourth  degree  polynomials  will  all  have  positive 
coefficients  and  will  factor  into  negative  real  and  ±  imaginary 
parts.  Thus  each  polynomial  will  appear  as: 


(s2  +  dj s  +  d2)  (s2  ■*  d3s  +  dj 


(81) 


The  graphical  technique  requires  solution  of  two  simultaneous 
equations:  (Reference  1) 


d 


i 


d  2  + 


C  d 
1  1  2 


(82) 


dt  = 


^  i  **/C  1 1  (d  2)  “  ^  1  2/C  1 1  (d  2  ) 

Cl5/Cn-  d22 


(83) 


Either  equation  is  rather  complicated  to  solve  by  plugging  in 
arbitrary  numbers.  For  stable  motion  it  is  known  that  values 
for  d.  must  be  positive,  d!  and  d3  will  be  much  smaller  than 
d2  and  d4.  The  latter  two  are  related  to  the  oscillating 
frequency  of  the  model.  The  above  type  of  characteristic 
equation  is  similar  to  the  linearized  equation  frequently  en¬ 
countered  in  studying  pitching  motion  of  low  speed  aircraft. 

In  these  equations  the  important  angular  quantities  are  pitch 
angle  and  angle  of  attack.  The  angle  of  attack  differs  from 
the  aircraft  pitch  angle  in  a  similar  manner  that  the  volute 
angle  of  attack  differs  from  the  cylinder  pitch  angle.  The 
mechanics  of  the  two  systems  are  different  since  the  volute 
angle  of  attack  is  the  result  of  aeroelasticity ,  whereas,  the 
aircraft  wing  angle  of  attack  is  purely  an  aerodynamic  phe¬ 
nomenon  (excluding  aeroelasticity) .  These  comparisons  mean 
that  the  volute  model  will  probably  have  a  short  period  mode 
and  a  long  period  mode,  as  do  aircraft.  The  frequencies  of  the 
two  modes  should  be  quite  different,  one  being  high,  the  other 
being  low,  analogous  to  an  aircraft's  short  period  and  phygoid 
mode  of  vibrations.  d2  was  chosen  as  the  parameter  easiest  to 
estimate  since  if  represents  a  frequency  quantity.  If  the  low 
frequency  vibration  is  estimated  to  be  about  2  cycles/second, 
d2  =  157.75.  This  is  a  good  approximation  to  begin  with  in 
solving  Equation  (82).  The  following  illustrates  the  technique 


Example : 

Solve 

for  the  roots  of 

Model  (3+3) 

(83) 

(82) 

d» 

d2 

di 

d2 

0.193 

160 

174 

165 

0.195 

169 

28.49 

168 

4.69 

168.5 

1.846 

168.56 

1.372 

168.57 

0.898 

168.58 

0.796 

168.59 

0.577 

168.595 

neg 

169.000 

The  right-hand  side  shows  that  increasing  d2  causes 
dx  to  drop  in  value  and  to  continue  dropping  until 
it  becomes  negative.  Inserting  168  into  Equation 
(83)  gives  0.193,  which  compares  to  28.49,  indica¬ 
ting  that  the  solution  is  still  off  a  good  bit. 

At  169,  d i  computes  to  be  0.195,  a  minor  change 
from  168;  however,  the  right  columns  show  dx  pass¬ 
ing  through  zero  which  would  indicate  that  the  solu¬ 
tion  must  lie  between  168  and  169.  Both  of  the  above 
sets  of  columns  plot  out  as  straight  lines  and  no 
further  points  are  needed.  Figure  22  is  a  plot  of 
the  above,  the  intersection  being  the  value  of  dx  and 
d2.  Direct  division  yields  d3  and  d4.  The  factors 
are : 


d  =0.194  d  =0.095 

1  3 

d„  =  168.582  d  =  8,161.4 

2  4  ’ 


With  these  values  for  the  factors,  the  quadratic 
roots  follow  directly. 
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Sj  =  -0.097  +  2 . 0  6  7  i  ,  i  =  /-l 

s2  =  -0.097  -  2 . 067 i 
s .  =  -0.048  +  14 . 41i 
s„  =  -0.048  -  14.411 

The  particular  solution  for  the  two  equations  of 
motion  is: 

0=  Aj e  •°-097t  (sin  2 . 067t  +  <^)  +  A2e  "0,048t  (sin  14 . 41t  +  ih2) 
6=  A3e  "0»097t  (sin  2 . 0 6 7 1  +  ^3)  +  A4e  (sin  14 . 41t+ip4 ) 


At  first  glance  it  would  appear  that  there  are  eight  arbi¬ 
trary  constants.  There  actually  are  only  four,  however,  since 
the  Ai  and  tyi  are  functions  of  four  independent  constants  found 
from  the  initial  conditions  0,  0  and  6,6.  The  technique  for 
finding  the  Ai  and  i^i  terms  whereby  the  complete  expressions 
can  be  written  is  outlined  in  Reference  2.  Since  stability  is 
the  prime  concern  of  this  program,  the  technique  is  not  included 
in  this  report.  The  equation  is  composed  of  two  angular  quan¬ 
tities  damped  at  different  rates.  The  mass  of  the  volute  dras¬ 
tically  changes  the  motion  of  the  forebody.  This  is  manifest 
through  the  additional  high  frequency  term  in  the  0  equation. 

If  the  mass  of  the  volute  was  sufficiently  small,  this  term 
would  have  negligible  effect  on  0  .  The  motion  would  be  deter¬ 
mined  primarily  by  the  leading  term.  At  t  =  10  seconds  the 
long  period  mode  has  damped  to  38  percent  of  its  initial  value 
and  the  short  period  mode  has  damped  to  62  percent  of  its  ini¬ 
tial  value.  Wind  tunnel  testing  showed  that  the  model  did  in¬ 
deed  have  the  kind  of  motion  indicated  by  the  solution  equation. 
The  short  period  mode  was  easiest  to  follow  with  the  eye  and 
damped  out  first.  The  remaining  high  speed  motion  was  then 
evident.  Unfortunately,  this  high  speed  vibration  only  damped 
out  partially  and  continued  to  vibrate.  This  limit  cycle  can¬ 
not  be  predicted  by  the  simple  linear  theory.  At  this  point 
in  the  two-year  effort  it  is  safe  to  say  that  the  residual 
motion  results  from  interaction  among  the  pitch,  yaw  and  stretch 
motion  of  the  spring.  It  was  somewhat  surprising  to  see  large 
amplitude  second  harmonics  appearing  in  the  bending  of  flexible 
helix  volute.  These  waves  in  the  spring  were  quite  evident. 

The  spring  motion  was  unsteady  and  appeared  very  complex.  An 
approximate  analytical  method  will  have  to  include  two  dimen¬ 
sional  oscillations  (coning)  and  possible  second  harmonic  ef¬ 
fects  (three  degrees  of  angular  freedom).  This  type  of  analysis 
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is  beyond  the  scope  of  the  first  year  study  but  may  be  the  sub¬ 
ject  of  the  second  year  study. 

The  remaining  roots  of  the  flexible  models  were  solved  for 
the  method  of  the  preceding  example.  The  roots  of  these  models 
are : 


Model  (3+2)' 


s2 

s3 

s4 

Model  (3+3)’ 

Si 

s2 

s3 

S4 

Model  (3+5)' 


Model  (3  +  6)  ' 

Sl 

s2 

s  3 

s4 


-0.104  +  2 . 164i 

-0.104  -  2 .1641 
-0.043  +  25 . 549i 
-0.043  -  25 . 549i 


-0.109  +  2 . 201i 
-0.109  -  2.2011 
-0.048  +  17 . 9  9  8 i 
-0.048  -  17 . 9 98 i 


-0.115  +  2. 040i 
-0.115  -  2 . 040i 
-0.053  +  12 . 380i 
-0.053  -  12 .3801 


-0.095  +  1 . 550i 
-0.095  -  1 . 5 5 0 i 
-0.044  +  8 . 660i 
-0.044  -  8 . 660i 


With  this  analytical  information  it  is  possible  to  construct  a 
curve  relating  the  damping  power  of  the  tail  along  with  its 
stability  power  to  the  volute  length  ly.  (figure  23  through  26) 
These  can  be  used  with  caution  to  predict  the  frequency  and 
damping  characteristics  of  volute  cylinder  combinations  dif¬ 
ferent  from  the  samples.  Caution  must  be  exercised  because  of 
the  vast  number  of  .assumptions  used  to  derive  the  analytical 
expressions.  Also,  the  forebody  in  a  particular  problem  must 
have  about  the  same  AB  as  the  plotted  values,  which  was  about 
0.440.  It  is  best  to  refer  back  to  the  theory  and  begin  with 
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Figure  25.  lv  Versus  Frequency 
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the  basic  expressions  for  the  a^'s  and  derive  the  stability 
margins  for  each  case.  Even  when  this  is  done,  the  answers 
represent  only  a  first  indication  of  the  model  stability  cha¬ 
racteristics  and  should  be  interpreted  as  necessary  but  not 
sufficient  conditions  for  dynamic  stability. 
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SECTION  IV 
FINAL  DESIGN 


The  preceding  sections  served  to  outline  two  methods  for 
finding  the  static  and  dynamic  stability  of  volute  equipped 
cylinders.  The  first  method  for  rigid  volutes  only  showed 
which  of  the  volute  parameters  entered  into  the  stability  cri¬ 
teria  and  provided  a  base  for  the  more  complex  flexible  analy¬ 
sis  which  followed.  Because  of  the  anticipated  complexity  of 
the  flexible  analysis,  the  first  approach  there  assumed  a 
small  mass  volute  but  allowed  bending.  This  method  took  into 
account  bending  by  altering  the  rigid  expression  for  Cma  and 

(Cm  +  Cm*)  by  a  constant  factor  dependent  upon  the  geometrical 
and^structural  constraints  of  the  volute  tail.  The  final  analy¬ 
sis  examined  the  dynamic  stability  of  several  models  through 
characteristic  equations  extracted  from  the  two  angular  equa¬ 
tions  of  motion.  It  remains  now  to  tie  the  preceding  work  to¬ 
gether  so  that  a  systematic  method  can  be  used  for  a  given  pro¬ 
blem.  One  more  piece  of  information  must  be  provided  along 
with  the  aerodynamic  analysis  to  allow  this. 

All  aerodynamic  work  derived  points  to  a  best  design 
length  volute  for  a  given  cylinder;  however,  it  may  not  be 
feasible  to  fit  the  optimum  length  into  the  cylinder  cavity 
(Section  I) .  The  maximum  length  of  the  volute  here  is  deter¬ 
mined  by  the  spring  thickness  and  the  cavity  diameter.  In  some 
cases  the  design  length  obtained  in  the  J  method  may  be  larger 
than  the  maximum  length  capable  of  being  stored  in  the  cylinder. 
To  make  matters  more  complicated,  the  spring  constant  kv  ap¬ 
pears  to  be  inversly  proportional  to  the  volute  length  for  a 
given  cross  sectional  inertia  which  changes  the  optimum  lv. 

This  latter  fact  represents  one  of  the  most  difficult  para¬ 
meters  to  estimate.  That  is,  for  a  given  material  wound  in  a 
spiral  to  a  length  lv,  what  is  kv? 

During  the  study  program  it  was  discovered  that  kv  varies 
inversely  with  lv.  Only  a  limited  number  of  volutes  were 
tested,  but  these  followed  Equation  (84): 

kv  =  (Constant)/  1^  (84) 

The  spring  used  on  Model  (3+4)  '  was  taken  as  the  center  point 
to  evaluate  (Constant).  If  z  =  y  the  proportionality  con¬ 
stant  was  found  to  be  0.501.  If  a  given  spring  is  close  to  the 
mentioned  spring  (Appendix  II)  in  cross  section  material  and 
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length,  Equation  (84)  will  give  good  answers.  Large  excursions 
must  be  dealt  with  separately. 

For  a  volute  to  have  maximum  efficiency,  it  must  be  wound 
as  lightly  as  possible.  That  is,  viewed  from  the  side,  there 
should  be  no  space  between  successive  coils.  To  determine 
whether  such  a  spring  will  fit  into  a  given  cylinder  cavity, 
Equation  (85)  can  be  used. 


ntj  <  R 


(85) 


To  find  the  length  of  uncoiled  wire  needed  to  construct  a 
volute  with  given  lv,  use  Equation  (86): 


lv“  2tt  (R  -  (n  -  1)  t,) 

The  volute  length  ly  is  approximately, 


lv  =  nt 2 

Spring  steel  is  recommended  as  the  volute  structural  material; 
specifying  the  cross  section  would  be  at  the  discretion  of  the 
designer.  The  key  flexible  volute  [(3  +4)']  test  model  mat¬ 
erial,  cross  section  was  0.020  x  0.200.  For  this  material, 
Equation  (84)  can  be  used  to  find  kv.  Using  the  J  method  to 
maximize  Cma  with  respect  to  volute  length  lv,  variable  kv  = 
f(lv)  should  be  used  if  Equation  (84)  can  be  used  or  an  alter¬ 
nate  derived,  based  on  different  material  than  0.020  x  0.200 
spring  steel.  Variable  kv  has  a  profound  effect  upon  Cma- 
Cma  falls  off  with  lv  faster  for  kv  =  f(lv)  than  kv  =  constant. 
This  fact  adjusts  the  'design  lv  to  a  smaller  value.  Example  6 
was  carried  out  for  kv  =  constant.  If  Equation  (84)  is  used 
instead,  Table  XI  results: 


(86) 


(87) 
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TABLE  XI.  Cjnaf/Rad. 


,  *  V  0.501 

Let  kv  =  — ■== 

■Jh, 

lv 

kv 

.. 

<J> 

J 

Cmaf 

0 

0.5 

0.708 

0.3343 

0.0021 

-0.03485 

pph 

0.4090 

0.3237 

0.0337 

-0 .11798 

2.0 

0.3543 

0.3145 

0.0612 

-0 .15786 

0.3168 

0.3019 

0  .0988 

-0.19142 

3.0 

0.2892 

0 .2876 

0.1415 

-C .21744 

3.5 

0.2678 

0.2725 

0.1866 

-0.23600 

4.0 

0.2505 

0.2572 

0.2322 

-0 .24695 

4.5 

0.2362 

0.2423 

0  .2767 

-0 .25180 

5.0 

0.2241 

0.2281 

0.3191 

-0.25162 

5.5 

0.2136 

0.2148 

0 .3588 

-0.24740 

; 

6.0 

0.2045 

0.2023 

0.3959 

-0 .24000 

6.5 

0.1965 

0.1906 

0  .4310 

-0 .22929 

7.0 

0.1890 

0.1802 

0  .4621 

-0 .21855 

7.5 

0.1829 

0.1705 

0  .4910 

-0 .20599 

OO 

o 

0.1771 

0.1614 

0.5182 

-0 .18911 

8.5 

0 .1718 

0.1533 

0  .  5424 

-0  .17919 
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These  values  were  plotted  in  Figure  27  and  compared  to  the  re¬ 
sults  of  Example  6.  If  kv  cannot  be  found  as  a  function  of  lv, 
use  the  best  constant  value  (from  elastic  testing  if  possible). 

It  is  likely  that  the  first  step  in  designing  a  volute 
stabilizer  will  be  to  insure  rigid  static  stability.  The 
approximate  dimensions  of  the  best  design  volute  are  obtain¬ 
able  from  the  simple  rigid  analysis  if  the  spring  constant 
is  very  large  (  2  pounds/radian),  or  if  the  volute  is  short. 
Generally,  the  volute  stabilizer  will  not  have  ky  as  large  as 
2  pounds/radian,  probably  closer  to  0.300  or  0.400  pounds/ 
radian.  In  these  cases,  the  J  method  should  be  used  to  follow 
up  the  first  order  approximate  rigid  analysis  and  the  tail 
length  altered  according  to  the  results.  Dynamic  stability  can 
be  checked  for  this  configuration  best  in  a  wind  tunnel.  How¬ 
ever,  the  two-dimensional  angular  freedom  analysis  will  provide 
at  least  a  first  approximation  to  what  the  dynamic  stability 
characteristics  will  be  like.  Since  the  observed  limit  cycle 
motion  referred  to  in  the  text  appeared  to  result  from  pitch- 
yaw  coupling,  a  more  precise  analysis  will  have  to  wait  until 
completion  of  the  second  year  study. 
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SECTION  V 


CONCLUSIONS 


The  preceding  analysis , together  with  data  collected  at 
the  Eglin  Air  Force  Base  low  speed  wind  tunnel  and  free  flight 
data  collected  at  the  Von  Karman  Gas  Dynamic  Facility,  proved 
that  the  volute  stabilizer  is  a  feasible  and  practical  approach 
to  providing  a  compact  stabilizer  for  cylindrical  forebodies. 

The  rigid  and  flexible  body  analysis  yielded  a  convenient 
method  of  estimating  the  effects  on  both  static  and  dynamic 
stability  when  the  various  geometrical  constraints  of  the  volute 
shape  are  altered.  For  each  case,  the  analysis  provides  the 
necessary  and  sufficient  conditions  to  insure  static  stability 
and  thereby  proved  that  properly  designed  and  matched  volute 
tails  insure  absolute  static  stability. 

The  extensive  amount  of  wind  tunnel  testing  conducted  at 
the  two  facilities  clearly  indicated  that  all  test  volutes  con¬ 
structed  of  coiled  spring  steel  operating  at  M  =  0.2  to  0.6 
and  Rn  =  1  x  105to  1  x  106  suffer  a  certain  amount  of  residual 
vibrations  which  apparently  never  damp  out.  This,  as  shown  in 
the  analysis,  can  be  offset  through  proper  volute  design;  how¬ 
ever,  it  appears  that  these  limit  cycles  can  never  be  removed 
entirely.  The  effect  on  stability  can  be  reduced  to  a  point 
where  the  forebody  wanders  through  ±  5  degrees  arc.  The  ampli¬ 
tude  of  the  limit  cycle  motion  is  sensitive  to  the  lateral 
spring  constant  of  a  given  spring  tail  and  under  some  condi¬ 
tions  will  produce  (Cm  +  CmO  >  0  causing  divergence  up  to 
a  certain  maximum  amplitude ,  the  largest  observed  to  be  about 
±  30  degrees. 

Under  nearly  steady  state  conditions,  the  cylinder  equip¬ 
ped  with  a  volute  stabilizer  will  trim  out  at  a  trim  angle 
which  is  controlled  by  the  volute  length  and  the  eg  location. 
This  trim  angle  can  be  reduced  to  nearly  zero  if  the  cylinder- 
volute  combination  is  carefully  matched.  If  this  is  not  done, 
the  combination  can  be  expected  to  trim  out  at  an  angle  dif¬ 
ferent  from  0. 
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(The  reverse  of  this  page  is  blank) 


APPENDIX  I 


ANALYSIS 


The  1-d  rigid  volute  model  is  analyzed  in  detail  below. 
Because  this  particular  model  is  more  amenable  to  analytical 
treatment  than  the  remaining  flexible  models,  this  section 
identifies  and  treats  the  basic  volute  parameters  which  affect 
aerodynamic  damping  and  pitching  moment.  The  effects  of  eg 
movement  and  forebody  induced  tail  blanking  are  included 
separately. 

At  the  outset  of  the  analysis,  the  following  assumptions 
are  made: 

(1)  The  volute  is  attached  rigidly  to  the  forebody 

(2)  The  volute  is  a  complete  cone 

(3)  All  area  behind  the  eg  acts  as  a  tail,  all  area 
in  front  of  the  eg  acts  against  the  tail 

(4)  Normal  force  coefficients  are  based  upon  cone- 
cylinder  combinations  and  strictly  cones  or  cylinders 

Volute  and  forebody  geometry  will  be  identified  according 
to  the  following  diagram: 


(Aero)t 


x 


A  stabilizing  pitching  moment  is  set  up  primarily  because 
of  a  normal  component  of  velocity  relative  to  the  volute  fore¬ 
body  axis.  By  finding  this  component  of  velocity  and  an  appro¬ 
priate  normal  force  coefficient  evaluated  at  ir/ 2 ,  the  normal 
force  can  be  calculated. 

Relative  Velocity 


Moving  out  to  a  point  on  the  forebody,  the  relative  velocity 
there  can  be  found  with  the  aid  of  the  vector  diagram  below. 


V 


r 

b 


(Vot  -  exbsin0)2  +  (0xbcos8)2 


The  angle  between  the  free  stream  velocity  and  the  relative 
velocity  obtained  from  the  diagram  is: 


<t>b  =  tan 


9xbcos0 
Vw-  0xbsin9 


Tail  quantities  are  obtained  in  a  similar  manner.  The  rotation 
induced  field  adds  to  the  free  stream  field  to  give: 


Vr  = 


(Vo,  +0xtsin0)2  +  (0xtcos0) 


4>t=  tan 


0XbCOS0 


+  0xtsin0 


Assumptions : 

(5)  Approximate  xb  and  xt  with  the  average  values  xb  and 

(6)  V„  >>0x  so  that  $b  *  0xb/VM,  <frt  s  0xt/Voo 


To  find  the  normal  component  of  velocity,  the  angle  of  attack  a 
must  be  found.  This  is  done  with  the  following  diagram: 


<x 


With  the  angle  of  attack  for  both  the  "b"  and  "t"  sections 
known,  the  normal  component  of  velocity  is  expressed  as: 


,rn  ,,r  - 
Vb  =  Vbsmab 

,.n  x,T  ■ 

Vt  =  Vtsinat 

Vb  =  Vasins  -  0xb 
=  VcoSine  +  0xt 

Moment  Balance 

To  find  the  total  aerodynamic  moment  acting  on  the  total  body, 
moments  are  summed  about  the  center  of  gravity. 

Assumptions : 

(7)  Use  average  moment  arms  h^  and  h-t 

(8)  Use  Cnto  based  on  cone-cylinder  values  evaluated 
for  ir/  2  radians. 

(9;  Sin 3  »  0 

(10)  There  is  no  tail  forebody  interference 

Then , 

-10  =  TpCntoht(Vco0  +  xt0)jVro0  +  xt0|St 
TP^nb0^b  (^oo9  '  xb0)Koo0  '  xb9ISb 
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Assumption : 


/\  * 

(11)  The  ratio  |6  -  Jke | / | 6  +  Jle |  s  l 

V  00  V  00 


Then, 

-10  =  7pCnt^htXtStV2oo 


(-'nb0^b^bxb 

Cnt0^tStxt 


2P^nt0^tStV2oo 


e  |  e 


Xt  •  I 

+  -Le 

VM 


Cnb^bSb 

CntJtSt 


Let 


vh  ■  hbsb/htst  •  «  -  9  *  *t/vJ 


Assumption : 
(12) 


C  nt^txtsj 
o 


1  + 


Cnb0xb 

Cntn*t 


Vh 


0 


^nt0^txtst 


*"nboXb 

1  ♦  _l°^vh 

Cnt0xt 


0 


and 


Cn^btSt 


Lnb0 

Cnc0 


Vh 


0  |  = 


Cntn^tSt 


Cnb0 

1  -  p — ”Vh 

Lnt0 


0 
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(13)  Let  Cnt  |e|S  Cnto9.  Static  normal  force  coefficients 
can  be  expressed  as  Cnt  =  Cnta  a  +  Csin2a.  Thus  Cnto 
CntctCtf/z)  +  C.  This  assumption  averages  Cnt0le!  over  9> 
and  requires  that  Cn*  be  small* 

1  I".  ^tcnt0^t^t  [*  cnb0xb  d~g 

-  19  =  7pV2  dS- 9  - 5 -  1  +  - “vh  .,0 

2p  “c*  d2S  Cnt0xt  VM 

L  cl-  J 


Cn«ohtSt  [l  -  — vh|  §  • 
dcS  1_  Cnt0  J 


From  linear  aerodynamics: 


19  =  |pV200dS  9  „  (Cmq  +  Cmi)  +  Cma9 


By  identification  of  common  terms : 


2StCntohtxt  ,  .  cnboxb  - 

<Cmq  *  <W>* - d*S  h  9 

CL 


Cnb0xb 


stcnt0ht  ^  cnb0 


■Vh  e 


vh'  =  htSt/dcS 


(cmq*cmS)=-cntov,;^ltggvhje 
-  -  cntovh'  [i  -  ^vhJ  e 
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The  effects  of  eg  movement  can  best  be  determined  by 
relating  the  movement  to  pitching  and  damping  derivatives 
through  the  above  equations.  Taking  the  pitching  moment  slope 
coefficient  first: 


Cnb^bSb  g 

cnt0htst 


for , 


Cnb^^bSb 
_o -  <  l, 

cnt0htst 


Cma  <  0 


This  constraint  can  be  related  to  the  eg  location  by  expres¬ 
sing  the  individual  terms  as  functions  of  xcg .  When  this  is 
done :  & 


"cma  =  ^  g  ^T^nb0  CSf  -  dcxcg)  (lb  -  xcg)  + 

\ 


CnvQSv(lb  "  xcg  +  alv)  -  Cnb0— ^&(xcgdc) 


The  individual  terms  as  functions  of  x 


are : 


eg 


Cnb0  =  Cnbo  0 


cnt0st^t 


Sb  =  ^cxcg 

Cnto  =  §7  (cnb0sb  +  cnv§vj  9 
^nbo^b  (lb  xcg  )t  +  ^nyf>v(ib  "  xcg 


+  alv) 


Sb  =  sf  - 


dcxcg 


Note,  that  because  Cnt0  is  n’-ade  up  of  normal  force  coeffi¬ 
cients  of  both  the  cylinder  and  cone  portions  aft  of  the  rota¬ 
tional  axis,  it  is  a  function  of  eg  location. 

After  some  reduction: 

-cma  =  xcg  (sf(Cnbo  +  CnvJ  -  Cnv^J  + 

^nvj1  ‘  Sf)  ^alv  +  lb J  +  T^nb0^f^b 

Where 

X  eg  =  X  Cg  /dC  >  Sf  =  Sf/S,  ly  =  ly/dc  ,  If,  =  lb/ dC 

A  more  convenient  form  of  the  above  was  obtained  from  Reference 
3  wherein  a  method  is  derived  which  allows  prediction  of  Cma 
about  any  axis,  if  Ct%  is  known  about  any  other  axis  in  the 
forebody.  The  transfer  equation  is: 
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Cn  has  to  be  determined  for  each  volute-forebody  configura¬ 
tion.  Cm^  can  be  obtained  by  setting  xCg  io  zero  in  the  eg 
equation  or,  directly,  from  the  definition  of  Cma.  When  this 
is  done : 


The  damping  derivative  turns  out  to  be  a  quadratic  function  of 
the  eg  location.  The  equation  is  too  unwieldy  for  general 
usage  but  is  duplicated  below  for  completeness: 


‘q  +  Cm«  dcS'  x2  eg  (<*  ^nbo1*1  +  z1v  enVo/2  +  a  ^bo ) 

;.g.  ^  "IT  ^nbQ1b2  ‘  VCnVoalv2  ’  Mb  CnVo/  ‘  a  Cnboj 


4  Cnb01b  +  1b  1v/2  Cnv0/2  +  a  ^nbQ  +  ^nv1v2lb 
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The  derivation  follows  that  used  for  Cma  as  a  function  of  eg. 
ic t i s  defined  to  be: 


xt  =  (lb  -  xcg)/2  +  alv 

xb  =  xcg/^ 


Since  Cma  was  found,  it  is  possible  to  use  parallel  axis  trans¬ 
fer  equations  to  find  the  effect  of  eg  movement  of  the  damping 
derivative.  This  equation  is  extracted  from  Reference  3  and 
appears  as: 


It  is  much  easier  to  use  than  the  one  derived  directly  from 
the  definition.  If  the  primed  parameters  are  measured  from  the 
nose  of  the  forebody  (consistent  with  the  definition  of  xcg) 
and  Cnq  is  small,  the  damping  derivative  as  a  function  of  eg 
location  is: 


for, 


|  ana , 


Cm 


a 


hrCnt 

§ 

dc 


Lcg 


0  (C 


cm^ 


+ 


x  eg 


0 
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htcnto 

<*c 


S\ 


0 


The  eg  location  for  maximum  damping  will  occur  at: 


x 


eg 


kht  nto  0 

Cna 


In  the  derivation  of  the  above  formulas  ,  an  average  moment  arm 
and  average  value  for  the  rotary  velocity  field  was  assumed. 
The  average  value  was  introduced  with  the  constant  multiplied 
by  lv(i.e.,  alv) . "a"  was  used  both  in  the  expression  for  aver¬ 
age  moment  arm  and  average  rotary^velocity .  The  value  will  be 
taken  as  •§-.  The  expressions  fgr  xt,  ht ,  and  later  xv^  are 
complex  and  depend  upon  6  and  0.  The  following  analysis  is 
devoted  to  this  subject: 


/s 


V2d  S  =  V2  s 
J  J  ave 

Sb  +  sv 


Plugging  in  the  necessary  values  and  reducing, 


.lo 

(V2O3sin20  +  ZVooXsinee  +  0  2x  2)dcdx 

'  0 

+ 

•  +  ly 

(V2msin20  +  ZV^xsinOO  +  0  2x2) 
lo 


~  (x  -  10)  dx  = 
lv 
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CT  dv1v+  dcl0)  (V2coSin20  +  ZVooSinee  xt  +  0  2xt2) 


A  2 


xt  = 


3  '*■0  +  3  ^0^-V  +  1V10  +  7 

21 n  +  lv 


,0=0 


xt  " 


Ip2  *  1V10  *  3  ly2 
210  + 


0  =  0 


*vt 

xy^  is  found  from  the  following  identity: 


v 


V 


2  -V 


x  dx  =  V2  ave  S 


Inserting  the  necessary  values  and  reducing, 


|v2oosin20  +  ZVooSinOCx  +  lo)0  +  02(x2  +  2xl0  +  102)t^- x  dx 
I  lv 

=  7  dvlv (V2ooSin20  +  2Voosin00xVt  +  02xVt2) 


cvt 


xVt2=  lv)2  +  41v1o/3  +  1o2»9  =  0 


xV-t  =  ^0  +  3  >  9  =  9 


lit 


ht  is  found  with  the  following  identity: 


JJ 


V2xdS  =  V2ave  h'^Sy 
S,r 
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Inserting  the  necessary  values  and  reducing, 


flv 

(x  +  10) 

0 

1 

2 


(V2cosin20  +  2Vcosin00(x  +  10)  +  02  (x  +  10) 
dvlvCV2oosin20  +  2Voosin00xVt  +  02x2Vt)h^ 


dv  , 

—  x  dx 
1V 


lnl 


0AV 


+  2102lv  + 


+ 


0 


0 


^t  =  ^0  +  3  =  ® 

The  above  forms  are  gotten  by  carrying  out  the  integrations  and 
alternately  setting  0  then  0  equal  tq  zero.  Th’e  three  variables 
span  a  range  of  values  which  can  be  considered  to  be  the  upper 
and  lower  bounds  of  each  during  one  cycle  of  oscillation.  It 
was  found  that  using  xt  =  l0/2  +  |-lv  gives  reasonable  answers 
and  is  easier  to  use  than  the  derived  expressions  which  them¬ 
selves  are  only  approximate.  ^h£  and  xVt  give  good  answers  de¬ 
fined  as  h£  =  (10  +  4-lv)  and  xVt  =  (lo  +  ylv)  •  Note  that  in 
the  expression^for  hi  at  0  =  0,tht  is  inversely  proportional  to 
the  square  of  Xyf  A  plot  of  ht  shows  that  for  certain  values 
of  lv,  h£  actually  decreases  for  increase  in  lv  up  to  a  point 
and  then  increases. 


The  preceding  forms  do  not  consider  forebody  interference. 
Because  the  forebody  interrupts  the  air  flow  over  the  volute 
tail,  the  effectiveness  of  any  volute  is  lowered.  To  find  what 
percentage  of  the  tail  is  blanked  out  by  the  forebody,  the  fol¬ 
lowing  approximate  analysis  is  offered.  Note  that  here  a  trape¬ 
zoidal  area  for  the  volute  is  used  instead  of  the  typical  tri¬ 
angular  planform.  This  is  done  because  at  low  angles  of  attack 
this  additional  area  has  a  significant  effect  on  the  volute 
performance  if  dy  is  large  or  about  Jv5.  The  following  diagram 
is  useful: 


dv 
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tandQ  = 


(y  ‘  dvr) 

j 

X 


(y  - 


x  tandQ 


Sv  (two  trapezoidal  pieces) 


ax  tand^x 


s  1  = 

°v 


tan  d0(lv2  -  x'2  ) ,  x*  = 


dV  “  dy' 


2(tan0  +  tand  0)  2(tan0  +  tand0) 


Sy  =  tan  d  o 


_  2 
-  d 


■4  (tan0  +  tand^ 


V  2  +  tan6o1v2 


Sy'  _  ly  dv 


dy2  2dy  lv8  (tan0  +  d  v) 2 


21, 


Sv"  (shaded  portion)  =  dy  (lv  -  x'  ) 


S  (total  area)  *  Sv*  +  Sv" 


Overshadowing  the  previous  analysis  is  the  obvious  fact 
that  the  volute  is  not  rigid,  but  flexible.  This  fact  alters 
the  previous  expressions  for  damping  and  pitching  derivatives. 
To  make  the  analysis  as  tenable  as  possible,  the  first  flex¬ 
ible  volutes  considered  will  have  small  masses  and  moments  of 
inertia . 

Using  the  following  diagram,  basic  flexible  volute  for¬ 
mulas  can  be  derived. 
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The  stability  equation  for  the  flexible  combination  can  be 
written  as  was  the  rigid  form  if  the  volute  contribution  con¬ 
tains  the  pitch  angle  and  the  volute  angle  in  the  aerodynamic 
force  expression.  Symbolically  this  is: 


(aero)b 


+  (aero)b  +  (aero)v  hi 


The  inertia  torque  of  the  masses  can  be  shown  (Reference  4) 
to  be : 


where , 


16  +  Kv6 

I =  'W^v  +  2MvRvrv  +  Iv*Ib 


Kv  =  MvRvrv  +  *v 


Written  out  the  whole  expression  is: 


16  +  Kv5  ='TpV200Sdc'| 


cnb0sbhb  unb0 


Cnh^S&hb  CnVoSvht 


dcS 


dcS 


dcS 


6  4 


2Cnb0sb^bhb  2Cnb0^b*b*1b  +  2cnVosv^vtht 

.  dr2S  dr2S  dr2S 


.  dc 
6  _  * 

2Voo 


CnV( 


Svht 


dcS 


6  <i>  + 


( 2CnVoSvXvb^  .  dc 


\  dc2S 


2VC 


r,  4=  1 6  +  6 


average 


To  take  a  closer  look  at  the  tail  reaction,  a  free  body  diagram 
is  drawn  and  moments  summed  relative  to  a  nonrotating  but  line¬ 
arly  accelerating  observer.  To  account  for  this,  the  inertia 
term  is  written  according  to: 


IM 


external 


m 


R  x  (co*R)  dm)  + 

y 


Mv (Rc .m 
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Which  yields  in  terms  of  problem  variables: 
l'i  +  MvRvrv0  ,  K  =  9+  6 


If  the  volute  is  idealized  to  have  negligible  mass  and  moment 
of  inertia,  the  above  reduces  to: 

-  kv6  =  (aero) v  a  lv  *  T pV2coCnVo  $  (6  +6)  a  1VSV 
This  is  a  constraint  on  6  expressible  as: 


6 


6 


CnvQ  $  a  lySv 

kv  +  TPV2OTcnVo  *  a  1vsv 

je 


Note  that  the  damping  force  was  assumed  small  compared  to  CnvQV 
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The  inclusion  of  this  expression  into  the  basic  flexible  form 
equation  reduces  the  expression  to  the  familiar  rigid  form  with 
variables  8,  0,  8.  The  coefficients  of  8  and  0  respectively 
are : 


18  +  Kv6  =  7pV200Sdc 


(  cnb0sbhb  Cnb  Svh't' 


drS 


cnb0sbhb  +  JCnVosvht 


dcS 


dcS 


+  7Pv2cosda 


2cnb0sbXbhb  +  2CnboSbXbhb 


dc  2  S 


dc  2  S 


2Cnv0SvxVbhb  2JCnv  Svxv  h^ 


dr2S 


dr2S 


fid 

=_c 

2V„ 


For  J  =  1  the  volute  is  totally  ineffective  and  its  contribu¬ 
tion  in  both  damping  and  stability  are  nullified.  For  J  =  0 
the  equation  reduces  to  its  rigid  body  form.  The  effects  on 
the  values  of  (Cmq  +  Cm^jand  Cma  can  be  inferred  by  reducing 
the  effectiveness ' of  the  volute  by  (1-J)  and  including  this 
value  instead  of  the  full  value  into  the  rigid  body  equations. 
By  doing  this,  the  tail  effectiveness  will  be  reduced  and  the 
overall  stability  of  the  combination  will  be  reduced.  Simi¬ 
larly,  the  damping  effectiveness  of  the  volute  suffers  by  the 
same  factor.  Note  that  since  the  whole  tail  (as  defined  in 
the  rigid  analysis)  is  not  affected  by  the  volute  flexing, 
the  aft  portion  of  the  forebcdy  which  helps  in  stability  is  not 
affected.  In  terms  of  J,  the  new  flexible  body  damping  and 
stability  coefficients  are: 


flexible 


Cma  [ 

8  jrigid 


$ 


JcnVnhtsv 

dcS 


3> 
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$ 


2JCnVohtSvXv 


If  the  volute  is  not  light  or  has  negligible  inertia,  the 
derived  forms  will  not  be  accurate.  A  complete  dynamic  analy¬ 
sis  is  somewhat  frustrating  since  the  results  are  so  compli¬ 
cated  as  to  be  almost  useless  for  general  design  purposes  un¬ 
less  a  convenient  means  of  interpreting  the  equations  is  pro¬ 
vided.  Here,  an  attempt  will  be  made  to  present  the  com¬ 
plete  equations  so  that  useful  design  information  can  be  ob¬ 
tained  . 

The  analysis  begins  with  the  general  equation  derived 
for  the  forebody  volute  combination,  and  the  equation  derived 
for  just  the  volute  without  the  small  mass  assumption.  The 
two  equations  have  the  general  forms  shown  below: 


10  +  K..6  =  -  a  0-a  0-a  6  -  a  6 

c  13  U  15  16 

(Iv  +  MvRvrv) 9  +  =  -  a2  30  -  a240  -  a2g6  "  a26^ 


The  coefficients  are: 


•  1  3 


rpV2ooSdc 


2CnbQSbxbhb  2Cnk^SbXbh|/  2Cnv^Svxv^.h-t 


dr  2S 


dr  2S 


dc2S 


dc* 

2V„ 


— '  “  '  OO 


14  ■  7PV2roSdc 


-n^  ^b^b  ^n^,  ^b^b 


cn  Svht 


dcS 


dcS 


dCS  J 


1  5 


7pV2  00Sdc 


2CnVoSvXvht  \  dc$ 


2V„ 


dr  S 


3  1  6  2  ooSdc 


cnVosvht^ 
dcS  J 
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a  =  7PV2coSd£  (  2Cnv0Svxvtalvyc^ 


>  a,.  =  2 pV2  Sd 


dc2S 


2VC 


2  4 


co^^C 


cnVo^va-*-v' 
dcS 


i  /  2Cn  Svx.,aL,  \  dr$ 

a  =  -  v2  q H  nyo  v  v  v  c 

2  5  2  PV  oo^ac  \ - 5 -  )  — — »  “26  "  2pv  »JUC 


dr  S 


2V„ 


.  -  7oV^Sdrf’Cnv°SvalI  1  ♦  *  k. 

dcS 


The  two  equations  for  which  a  solution  is  sought  are 


a: i  +  ai26  +  a i a 0  +  ai46  +  ais6  +  ai66  =  0 


a  2 1 0  a226  a  2  3  9  a  2  4  0  3 2  5 5  a2  66  =  0 


To  find  a  solution  assume  that: 


9=0  est,  6=6  est 


Insertion  into  the  equations  of  motion  yields: 


S2(all0  +  a126)  +  S(a13§  +  als6)  +  (.a140  +  alg6)  =  0 

S2  Ca210  +  a22^  +  Sta2  30  +  -:25^  +  ^ 3  2  4 0  +  a26^  =  0 

These  equations  for  the  assumed  <5  and  0  will  give  solutions 
different  from  zero  only  if: 


fa  x iS2  +  a13S  +  a14)(a12S2  +  a15S  +  a16) 
(a2  iS  +  a  2  3S  +  a.  2  (.3-  %  2^  +  a  2  5  ^  +  a26^ 


This  requires  that: 


CnS'  +  C12S3  +  C13S2  +  C14S  +  C15  =  0 


=  [0]  for  M  e  M 


C1X  =  f 3  1  1  3  2  2  '  a!2a2l)>  C  1  5  =  f 3  1  4 3  2  6  '  a  l  ^  2  ^  etC 
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Solutions  to  the  characteristic  equation  are: 


II 

f* 

CO 

-n,  +  ip ! » 

(1 

04 

CO 

-nx  -  ip 

S3  " 

* n  2  +  ^2 

S4  " 

-n2  -  ip 

In  terms  of  trigonometric  functions: 

0  =  e‘ri1t(A  cos  p  t  +  A  sin  p  t)  +  e'n2t(A  cos  p  t  +  A  sin  p  t) 

1  1  2  1  3  r  2  4  v  i  1 


6  =  e"ni^-(B  cos  p  t  +  B  sin  p  t)  +  e"n2^(B  cos  p  t  +  B  sin  p  t) 

112  1  3  2  4  2 

Before  delving  deeper  into  the  solution  equation,  it  is 
useful  to  examine  the  characteristic  equation  further  since 
stability  of  the  system  can  be  found  by  examining  its  co¬ 
efficients  and  factors.  Static  stability  can  be  inferred 
immediately  by  examining  the  values  of  the  constant  portion 
(C i s )  of  the  characteristic  equation.  ft  should  be  positive. 

If  it  is  positive  for  one  configuration  and  a  single  design 
parameter  is  varied  so  that  C i s  changes  sign,  then  one  diver¬ 
gence  (static  instability)  appears  in  the  solution  (Reference 
1  ) .  To  find  out  what  effect  this  change  has  on  dynamic 
stability,  the  characteristic  equation  needs  to  be  analyzed 
further.  Routh's  Criteria  is  a  convenient  means  for  finding 
whether  the  system  is  dynamically  unstable.  This  method  is 
outlined  in  Reference  1.  Briefly,  the  technique  is  as  follows: 

The  quartic  is: 


Cn^  +  C12s3  +  C13s2  +  C14s  ♦  C15  =  0 
Form  two  rows  as  follows: 


r  r  r 

^ 1  1  ^ 1  3  5 


I  C  i  2  C !  4  0 
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Construct  additional  rows: 


P31  P  3  2  0 


P4!  0  0 


where : 


P  3  1  =  1  2  ^  1  3  '  ^  l  1  ^  1  4  »  P  3  2  ^12^1 


p  =  p  r  -pc 

1  4  1  3  1  1  4  3  2^12 


For  the  system  to  be  dynamically  stable,  each  element  in  the 
first  column  must  be  positive.  If  the  mentioned  change  in  a 
single  design  parameter  results  in  an  element  in  the  first 
column  changing  to  a  negative  value,  then  the  system  will  have 
one  divergent  oscillation. 

Once  the  volute-forebody  characteristic  equation  has  been 
obtained,  any  of  several  root  extraction  methods  can  be  applied 
to  determine  the  degree  of  damping  and  the  two  free  vibrations 
of  the  system.  It  is  intuitive  that  the  system  will  vibrate  in 
both  modes  and  will  be  lightly  damped  so  that  the  solution  roots 
will  have  negative  real  parts  and  imaginary  parts.  Since  the 
equation  is  fourth  degree,  there  must  be  two  quadratic  factors 
obtainable  by  the  following  method: 

Rewrite  the  quartic.as: 

(s2  +  dxs  +  d2)(s2  +  d3s  +  dj  =  0 

Expand  this  equation  and  equate  coefficients  to  the  original 
quartic : 


C  1  3 

d,  +  d  .  - — 

1  3  ’  f 

U  1  1 


d  +  d  d  +  d 

2  13  4 


d  3  +  d  2 , 


d  2  d  4 
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These  equations  can  be  solved  to  yield: 


d  =  12  ± 

u  1  - 


2Cm 


IS  2 


C i  id 2 


di  " 


Cl./Cll(d2)  ‘  C12/C11(<1!2) 


'  2 


A  plot  of  these  two  equations  will  reveal  the  common  points  at 
one  or  two  intersections.  These  are  then  the  factors  of  one 
quadratic,  the  other  two  following  directly.  The  equations  are 
analyzed  below: 


m/c  "  C  1 2 / C . 


di  - 


^11  >  ^1  2  >  ^13  ’  ^14’  C1  S  >  ^ 


c»/cl:  d2: 


d2  -  dj  -  C12/Cii 


di  -*■  00 ,  d. 


-+  C  t  5 


dl  <  0»d2  '"C14/C1 

Looking  at  the  second, 


d  i  = 


2C i  i  \  2C  j  i 


iilV  -  ^+d  Jill-  T 


Cll  2  C^d; 


±  -  - 

2 


ci2  V  c13  ,  ci5  "7 

- -  -  +  d,+  - 

•C  i  i  /  C1X  Cjjdj 


l15  -2 

1  -  —  d. 


v*  1 

c  2 
u  1 1 


Maximizing 


C  l  s  ;2 


C  =  1  *  -  d2  ,  d2  =  ± 

Cm 
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APPENDIX  II 


DATA 


The  data  portion  of  the  program  was  conducted  in  two  parts. 
The  first  consisted  of  free  flight  data  obtained  from  the  Von 
Karman  Gas  Dynamics  Facility.  Static  and  dynamic  data  was 
obtained  at  the  Eglin  Air  Force  Base  low  speed  facility  at 
Mach  0.19  for  Rn  =  5x10s.  The  free  flight  data  is  contained 
in  the  report  "Volute  Stabilization"  project  VA0086. 

Since  the  volute  is  a  coiled  flexible  spring,  it  was 
necessary  to  approximate  the  rigid  model  of  a  volute  with  a 
solid  body  grooved  in  a  staircase  fashion  to  represent  the 
volute  coils.  Five  aluminum  and  one  epoxy  models  were  construc¬ 
ted.  Three  cylinder  forebodies  were  made  to  mate  these  volutes. 
The  test  bodies  with  complete  geometrical  and  physical  proper¬ 
ties  are  shown  in  Figures  1 1 - 1  throueh  II-3.  These  models 
formed  the  core  of  the  rigid  body  data  bank.  Several  attempts 
were  made  to  construct  flexible  models  from  these  by  inserting 
a  leaf  spring  between  the  forebody  and  the  volute,  giving  a 
two-degree  angular  freedom  model.  Unfortunately,  the  models 
generally  failed  structurally  because  the  aluminum  volute  was 
too  heavy.  Some  flexible  data  was  obtained  at  Eglin  Air  Force 
Base  for  each  flexible  model.  The  coiled  spring  tails  on 
models  (1-2)  are  shown  in  Figure  1 1  -  4 . 

Static  data  for  the  rigid  models  is  plotted  in  Figures 
1 1  -  5  to  1 1 - 1 0 .  Flexible  data  for  flexible  models  is  shown  in 
Figures  1 1  - 1 1  and  11-12. 

Two  specialized  pieces  of  equipment,  an  accelerometer 
stand  and  a  dynamic  shaker  stand,  were  developed  for  the  test 
program  at  Eglin  to  obtain  dynamic  derivatives.  Neither  of 
these  items  were  fully  tested  during  the  program  because  of 
time  constraints  and  may  be  used  extensively  during  the  second 
year's  work.  The  operation  of  each  is  described  below: 

Accelerometer  Stand 

The  accelerometer  stand  was  developed  for  mounting  in 
the  Eglin  low  speed  tunnel.  Inside  the  strut  are  assor¬ 
ted  electronics  for  transferring  a  voltage  signal  from  an 
accelerometer  model  to  a  special  discriminator.  The  out¬ 
put  from  the  discriminator  drives  an  oscillograph,  giving 
a  trace  of  the  acceleration  of  the  volute  versus  time. 
Diagrams  of  the  strut  and  electric  package  are  provided 
in  Figures  11-13  and  11-14,  respectively.  With  the  various  filter  modi 
fications,  the  strut  may  be  used  in  flexible  test  work 
during  the  second  year's  work. 
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Dynamic  Stand 


The  dynamic  stand  enables  extraction  of  the  basic  aero¬ 
dynamic  derivatives  of  a  given  model  by  forcing  it  into 
resonance.  Data  is  obtained  with  a  high  speed  camera  and 
reduced  according  to  a  scheme  outlined  in  Reference  4  . 

The  stand  is  illustrated  is  Figures  11-15  and  11-16. 

A  sample  of  the  free  flight  data  obtained  from  the  Von 
Karman  Gas  Dynamic  Facility  is  included  as  Figures  11-17  through 
11-20.  Reduced  data  for  the  aerodynamic  coefficients  has  been 
extracted  from  "Volute  Stabilization"  report  and  included  for 
completeness  and  comparison  with  the  planar  oscilatory  data 
obtained  at  Eglin  Air  Force  Base. 

Models  are  formed  by  combining  the  cylinders  with  the 
volute  tails.  Their  designation  is  provided  by  identifying 
each  model  by  the  cylinder  forebody  number  and  the  volute  tail 
number,  such  as  3+2,  which  is  a  model  having  cylinder  3  as  its 
forebody  and  volute  2  as  the  tail.  Some  flexible  tests  were 
conducted  by  inserting  a  leaf  spring  between  the  cylinder  and 
volute.  These  models  are  identified  as  above  except  that  they 
are  primed  thus:  (3+2V .  The  physical  properties  used  to  ob¬ 
tain  the  basic  inertia  terms  are  shown  in  Table  1 1  - 1 . 


TABLE  II-l.  VOLUTE  MODEL  MEASUREMENTS 


Model 

rv 

Rv 

*v 

I 

3  +  1 

0.05533 

0.1516 

0.00003636 

0  .0002044 

3  +  2 

0.10033 

0.1516 

0.00011210 

0.0003428 

3  +  3 

0.14408 

0.1516 

0.00027150 

0.0006050 

3  +  5 

0.21592 

0.1516 

0.00073500 

0  .0015090 

3  +  6 

0.38850 

0.1516 

0.00422200 

0  .0057600 

The  reduced  data  from  the  Von  Karman  test  program  can  be 
used  for  design  purposes  as  well  as  providing  an  extensive  data 
bank.  Various  aerodynamic  coefficients  were  reduced  from  the 
raw  data  for  a  number  of  models,  distinct  from  the  ones  used  at 
Eglin.  Table  1 1  -  3  lists  the  data  from  the  Von  Karman  Gas 
Dynamic  Facility.  To  assist  in  correlating  that  data  with  pre¬ 
vious  model  configurations,  Table  1 1  -  2  is  provided. 
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VOLUTE  4 


1* — 0.257 


3.448 


Figure  1 1  -  2  .  Rigid  Volutes  :  4,  5,  and  6 


1.620 


1.820 


3.440 

Model  Pitch  center  locations  a 


Figure  1 1  - 3  . 


CYLINDER  1 


CYLINDER  2 


CYLINDER  3 

Weight  (Pounds) 
0.434 

IB  (Slugs  -Ft?) 
0.00009558 


the  same  as  the  cylinder 


0  20  40  60  80  100 

Pitch  Angle  9  (Degrees) 

Figure  1 1  -  5.  9  Verses  Cn  For  Rigid  Volutes 


103 


r 


Pitch  Angle  9  (Degrees) 

Figure  1 1  - 6.  9  Versus  For  Rigid  Volutes 
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0  0.05  Time  (Seconds)  0.15 

Figure  11-18.  Von  Karman  Data  Sample,  0.75  Inch  Cylinder  With  Small  Volute 


(seajtSaa)  g  ojguv  tpiTd' 
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0  0.05  0.10  0.15  0 ; 20 

Time  (Seconds) 

Figure  11-19.  Von  Karman  Data  Sample,  1.25  Inch  Cylinder  With  Small  Volute 


(soipui)  TOABJJ,  ITJOp^aSA  puv  TB1UOZIAOH 
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0  0.05  Time  (Seconds) 

Figure  11-20.  Von  Karman  Data  Sample,  1.25  Inch  Cylinder 


TABLE  1 1 -2.  CYLINDER  MODEL  CHARACTERISTICS 


Cylinder  Length 
(Inches) 

Weight 

(Pounds) 

Moment  of  Inertia 
(Slugs-Ft2) 

xcg 

i  (Inches) 

0.75 

0.059 

2.2  x  10"‘ 

0.45 

1.00 

0.079 

r~374  x  10'b 

0.55 

1.25 

0.088 

4.3  x  1 0  '  b 

0.72 

2.08 

0.380 

1.36  x  10-3 

1.25 

2.78 

0.500 

2.05  x  10‘3 

1.65 

3.48 

0.610 

2.95  x  lO'4 

1.96 

Table  I I - 2  lists  only  the  cylinder  forebody  without  the  volute  tail. 
Only  two  different  type  tails  were  used  for  the  Von  Karman  tests.  These 
two  types  are  shown  in  Figure  11-21.  The  smaller  of  the  two  was  used  on 
cylinder  l.agths  0.75,  1.00,  and  1.25  inches.  The  larger  volute  was  used 
on  the  remaining  larger  cylinders.  Both  volutes  were  constructed  of  epoxy 
or  silicon  rubber  corresponding  to  rigid  and  flexible  volutes,  in  both 
cases,  the  volute  tail  was  hollow  to  give  a  more  realistic  value  for  Iv. 

This  is  especially  important  in  order  to  have  meaningful  free  flight  data. 

During  the  tests  with  the  above  mentioned  volutes,  it  was  decided  to 
use  a  number  of  coiled  spring  tails  since  these  were  available  in  one  size 
(Figure  I I -4 , Flexible  Volute  2).  These  small  spring  tail  volutes  stabilized 
the  smaller  cylinder  forebodies  as  previously  mentioned.  Data  listed  in 
Table  I I -3  pertains  to  this  particular  tail  for  configurations  075-1,  100-1, 
and  125-1. 


Figure  11-21.  Volute  Tails  Used  In  Von  Karman  Test  Prograir. 
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TABLE  I I- 3.  VOLUTE  STABILIZATION  DATA  (Continued) 
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TABU;  11-3.  VOLUTE  STABILIZATION  DATA  (Concluded) 


rH 

CD  T) 

E  3 
CJ  OO 

1 

1 

1 

i 

i 

I 

CM 

1 

LO 

1 

1 

1 

• 

i 

i 

H- 

• 

*3- 

i 

\D 

to 

i 

i 

i 

i 

o 

• 

rH 

1 

o 

• 

rH 

o 

rH 

1 

-2.55 

-1.33 

-1.22 

• 

CC 

1 

— 

— 

E  Lh 

1 

i 

vO 

1 

l 

t 

i 

i 

i 

1 

l 

o 

at 

1 

1 

cru  T2 

.» 

i 

1 

I 

r 

i 

i 

i 

1 

l 

rH 

1 

1 

s  « 

u  +  c- 

1 

I 

1 

i 

t 

i 

i 

1 

l 

1 

1 

1 

1 

Q 

\D 

co 

at 

o 

CM 

CM 

i 

vO 

00 

m 

oo 

at 

LO 

U 

• 

• 

• 

• 

• 

• 

• 

« 

i 

• 

• 

• 

• 

• 

• 

LO 

CM 

rH 

co 

LO 

CM 

CM 

i 

CM 

to 

CO 

CM 

CM 

CM 

N 

O  V) 

JO  3  4; 

U  4)  P  <U 

O  H  (-1  P 

m 

LO 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

m 

o 

o 

3  oo<  oo 

rt 

oo 

00 

00 

oo 

at 

at 

3  3  4) 

rH 

rH 

rH 

rH 

-J  <  iwQ 

Ow 

O 

rH  i 

oo 

\D 

o 

to 

to 

tO 

rH 

rH 

00 

CO 

nj- 

co 

o 

m 

1  o 

3 

Q 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

00 

co 

8  •  rH 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

t 

• 

• 

4i  P 

H 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

BO  to  X! 

o 

Ol 

Cl 

at 

at 

CT) 

00 

00 

at 

t- 

r-M 

•o 

in 

LO 

LO 

LO 

in 

LO 

in 

to 

CO 

to 

co 

co 

CO 

co 

s 

• 

• 

♦ 

♦ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

» 

t 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

1 

3 

oo  3 

rH 

rH 

o 

o 

rH 

rH 

rH 

rH 

rH 

rH 

o 

o 

o 

rH 

rH 

•H  O 

1 

1 

1 

1 

l 

1 

i 

1 

l 

1 

1 

1 

1 

l 

l 

<H  *H 

o 

o 

in 

o 

o 

o 

O 

o 

o 

O 

o 

o 

oo 

o 

o 

CO  P 

o 

o 

r- 

o 

O 

O 

o 

o 

o 

O 

o 

o 

o 

o 

O  3 

rH 

rH 

o 

rH 

iH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

CM 

rH 

rH 

U  P 

o- 

3 

o> 

o 

iH 

CM 

CO 

to 

rf 

LO 

tO 

r^. 

oo 

at 

o 

rH 

O 

LO 

LO 

LO 

in 

in 

in 

in 

in 

in 

LO 

in 

LO 

to 

to 

P 

u 

t/> 

P 

•H  O 
00  •> 

•H  *  O 
*3  DO 
H  • 

4)  Jl  r-l 
<U  .H 
P  X  « 
j:  4>  m 
HHN 
4-1  • 

P  O  • 
(/>  rH  <U 

P  t/l  P 
•H  (/>  p  3 
4n  <U  4>  rH 
P'3  o 
4>  3  CS  > 
rC  O  -H 
P  *rl  r-H  p 
•3  >s  <U 
P  G  O  bO 
<U  *H  P 

p  g  3 

n  n  i~i  o  rH 

p  t/l 
4)  3  *3 
P3  D 
G  to 
•rt  <u  3 

rC 

-P  4) 

4>  p 
rH  00  4) 

J3  G  £ 

3  *H 
P  JH 
O  CM 
4)  rH  I 
>HH 
O  O  t-H 
JO  4h 
3  4) 

P  P 
4)  -rl  3 
JO  00  00 
P  -H  -H 
tJ  U, 

C 

•rt  4>  G 

JO  *rl 
4)  P 


N  4) 

•H  *4-> 

V)  JO  3 

P  rH 
H  MO 
4)  (0  > 
T>  <D 

ohh 

S  rH 


P  3 
<u  4)  g 
JO  T)  Ml 
p  CO 
•H  4) 
TO  rH  jo 

CO  OnH 

•H  O 
4h 


CM 


<D  •  • 

'  O  JO  no  rH 
H  P  -rH 
-  OOT) 
W‘H  0 
3  p  3 


122 


5  in  length.  The  larger  cylinder's  used  the 
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